— 





Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1985 


An analysis of symmetric reinforcement of 
graphite/epoxy honeycomb sandwich panels 
with a circular cutout under uniaxial 
compressive loading. 


Sullivan, Patrick D. 


http://hdl.handle.net/10945/21158 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| г D U DLEY research materials and institutional publications created by the NPS community. 
FW узу, Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


WW KN OX appointed — and published — scholarly author. 

OM LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 73 





http://www.nps.edu/library 


کیہ ہہ ا او سیا میسو Кдра‏ 
а iy‏ رو را چرام 
ماس سے سس مو ےید یر چیہ 
tcov icr mr M e i-o eg tmt‏ 


Po aids dpi UA оф ВЕРНИ ОТ радар арфа A 
ер ды سار رس نین‎ a ا‎ ара پچ تیر سید سای‎ etude ا ا ا ا‎ 
اھ جوا مان ابو‎ O و‎ ндеу A, REN ач аен ар 
4 Ku ће) че иа 
مخ‎ apii را‎ зе رر انور‎ PRA OO cae Or 





- اس 
"v‏ مد وہ В ah an o.‏ او ووی ییون ا ا К Сс Баа Gp sn‏ ا ا ا m hamatum sumi‏ سا سا و ای 
e‏ ید У борави н пора о Врања‏ دہ مر دو حر ss‏ 
АЕ њи — а a қ е‏ رھ ی چ وک е C‏ 
APP AGLI EPIL abet ој „л npg ұғыл алу бы ча ч ٠ 8‏ 
правио аи d bad qm qti ap do dio te‏ بت فی ضرم р dachte‏ 
E күлер оов ~ e‏ وھ و کی ت رہ ی ب o‏ ا = ر ی ب 
IN 7 ду roo qp ph nm em nua sei enis доње На d t dde v pads‏ ~ 
a a A aa H +‏ اب س وا di наи‏ وو وچ چچ چب ерем» тағынған ата tin ete‏ 
ро ~ u дікі = eb.‏ 


nd 


m M 
Py atte od Рака di ا سی‎ 25 x " = v 
ld куы ы" = а 71 d T i T 3 Ре dumm f s "m 7а e + 
РРА И m aua P apo ih go Y a. PPP DP n dili d dmn <a ФН rated terre ner part с 

а I - و یوک و ا‎ Т GT I ж УРДУ РРА а poer ч у 

3 ا‎ os ү рше сла өл aA И xwe o wee = pa End 
е” СР ДА Де d epp crm Perf poer ya y Povery ~ep er NG OO L i. d = 
ый 8 Š L... ie w^ ИРА оо اید‎ дај dad بد‎ и 

ма یم‎ k. n РА сан d 
ne y ғ. б: e v a < сре ври күчө ч е" Еш nias ۵ / fi қ 
N (n F До — ићи 

po: ы ~ ~ Тр Дей ри об 
Las 


” ” = и E 
` M ~ Жетес аға тары РР Ири и 

"ар nud 5۶ PP AE m€— к “if e - и и we enu ИК ла дады = s ۹ е 

дай ље = ص‎ A LE ore - ды سز‎ ч» Ры, ue Pam am mE aur e Pr AR ы Pto MEN КАНЕ مض‎ m ~ ара“ >" | У n" vgl ^g? e qf жай ы ЖА” na РА ыды 9. < ibs B = а 214 
E e eem ç <. um ree e Мс» Á Thi... PEEL" ME i iie ро a o vewe Pe Sa ap, ұлға - 3 
ہے‎ na”. p. ы К ШЕ „се кыбы [мич үү عو‎ E - ےراس‎ Б E s i А ی ی و‎ la iq овна Дн АЛЫНЫА ИН 
ымы жен - فو اا‎ metano tI, go aw POLE OOO e ` 4 3 м F 
у = „ “РХ ту 


=> - 
и ДОРА СЕКСА aaah Pia мех = 
babe x 


3 Ка Да < 

LI ... и - ~ Š чи ا‎ : 

Ë в : b $ 4 ےت ہا موم‎ ~ ум у m y w. —9À C АД. 

L Z. нв 3 6 -—- ñ E e^ cc te qt, sf Ig a rta ne ова ird Ф ьа >. S‏ وم سی یہ 

ыыы a МАРИА n a pq D manne چ‎ > aa. i. “Q... ae EA - - т A ЧК А2 4 — а سا‎ ot АД» УЕ ЖШС ا یا سمخ مہ بی‎ Ku... 

E ~ Док a ana iu, kasih. V ⁄ P^ ےی‎ ww К 3 b ін j s Ме مو‎ р kd 

~ wa w. `... 2 e КРЕ УЕА. ДАУ 4 ў КА аф” E : қ A 

В Е اص دوپ‎ + md rm - i 
P P e“ سے‎ e Parr aa Ie Дон ا‎ ~ 

=. = рече" 


е оно от ет 











NAVAL POSTGRADUATE SCHOOL 


Monterey, Galifornia 





AN ANALYSIS OF SYMMETRIC REINFORCEMENT OF 
GRAPHITE/EPOXY HONEYCOMB SANDWICH PANELS WITH 
A CIRCULAR CUTOUT UNDER UNIAXIAL 
COMPRESSIVE LOADING 
by 


Ра 8 u Sullivan 


December 1985 


Thesis Advisor: 





рр = for pubdwaewreleoase; distribution unlimited. 


hee gy | 













"URITY CLASSIFICATION OF THIS PAGE 
REPORT DOCUMENTATION PAGE 


REPORT SECURITY CLASSIFICATION 1b. RESTRICTIVE MARKINGS 


0 "аа. 


3 OISTRIBUTION/ AVAILABILITY OF REPORT 
Approved for public release; 
٠1-31676071100 


SECURITY CLASSIFICATION AUTHORITY 








DECLASSIFICATION / OOWNGRADING SCHEDULE 


го Mali ne خر اق‎ 


PERFORMING ORGANIZATION REPORT NUMBER(S) S. MONITORING ORGANIZATION REPORT NUMBER(S) 





NAME OF PERFORMING ORGANIZATION 6b. OFFICE SYMBOL Î 7a. NAME OF MONITORING ORGANIZATION 
(If applicable) 
E Postgraduate School Code 67 Naval Postgraduate School 


ADORESS (City, State, and ZIP Code) 7b. ADDRESS (City, State, and ZIP Code) 


Enterey, California 93943-5100 Monterey, California 93943-5100 


| 8b. OFFICE SYMBOL 
(If applicable) 


NAME OF FUNDING / SPONSORING 


9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION | 







| کے ے کیا کح اا‎ Ae ale islet WS ABUL lial à а 


ADDRESS (City, State, and ZIP Code) | 10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT Î TASK IWORX UNIT 
ELEMENT NO. NO NO {ACCESSION NO 


, de wh. 





Mi ee іы Ln 


TITLE (include Security ED 
N ANALYSIS OF SYMMETRIC REINFORCEMENT OF GRAPHITE/EPOXY HONEYCOMB SANDWICH 


ЕМЕЕ5 WITH A CIRCULAR CUTOUT UNDER UNIAXIAL COMPRESSIVE LOADING 


| PERSONAL AUTHOR(S) 
M Evan, Patrick D. i 


TYPE OF REPORT 13b TIME COVEREO [14 DATE OF REPORT (Year, Month, Day) |15 PAGE COUNT 
FROM TO 1985 December | 


по“ '= thesis 343 4 
SUPPLEMENTARY NOTATION | 
| | 

š 















COSATI CODES 18. SUBJECT TERMS (Continue on reverse if necessary and identify by block number) š 
FIELD Composite: reinforced composite; circular hole; i 
Ш 1 _} Graphite epoxy panel; compression | 
ABSTRACT (Continue on reverse if necessary and identify by block number) ч 
An experimental and computational analysis was made of stress/strain 
E ntrations around a reinforced ciruclar 1.00 inch diameter circular 
1tout in HMF330C/34 (cloth) graphite/epoxy (G/Ep) and fiberglass/phenolic ! 
E cemb sandwich panels under uniaxial compressive loading. The test 
pecimens were 10.00" x 8.50", eight ply quasi-isotropic ([0,445,90,core] ) 
ES Ihe reinforcement consisted of either one or two additional G/Ep 
Eres co-cured to the outside of each facesheet. Three general reinforce- 
ent configurations were considered: round, square and strips parallel to 
ne applied load. The analytical results demonstrated that small amounts 
- reinforcement could greatly increase the strength-to-weight ratio. The 
idication was that concentrating the reinforcement close to the cutout 
telded the greatest decrease in stress concentration. A program of ex- 
*rimental validation of the analytical results experienced some problems 


` سینا تضم سیا م wa. o e i o‏ فوبےر سے اص ساس 


Bal cts dada noti 


DISTRIBUTION / AVAILABILITY OF ABSTRACT 21. ABSTRACT SECURITY CLASSIFICATION 
CXUNCLASSIFIED/UNLIMITED ÜC] SAME AS RPT Clotic users $ Unclassified 

L NAME OF RESPONSIBLE INDIVIDUAL 22b TELEPHONE (Include Area Code) | 22c. OFFICE SYMBOL 
IU. Bank (408) 646-2582 Code 035451 


FORM 1473,84 MAR 83 APR edition may be used untıl exhausted SECURITY CLASSIFICATION OF THIS PAGE 
Allother editions are obsolete 


J: 









Dell ha‏ فة نن Д Анас:‏ ف همی 








SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered 


Block 3. (cone a 


with prematrue panel failure caused by the facesheets separating 
from the core. It generally confirmed the analytical rasu g 
however. Further experimental tests on promising reinforcement 
configurations are justified based on these results. Properly 
designed reinforcement around cutouts in composite panels can 
significantly reduce the stress concentration and holds the 
promise of far lighter and stronger aerospace structures. 


2 


SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 








Approved for public release; distribution unlimited. 


An Analysis of Symmetric Reinforcement of 
Graphite/Epoxy Honeycomb Sandwich Panels with a 
Circular Cutout Under Uniaxial Compressive Loading 


by 
Patrick D. Sullivan 


Commander, United States Navy 
B.S., United States Naval Academy 1969 


Submitted in partial fulfillment of the 
requirements for the degree of 


AERONAUTICAL ENGINEER 
from the 


NAVAL POSTGRADUATE SCHOOL 
December 1985 


ABSTRACT 


An experimental and computational analysis was made of stress/ 
strain concentrations around a reinforced circular 1.00 inch 
diameter circular cutout in HMF330C/34 (cloth) graphite/epoxy 
(G/Ep) and fiberglass/phenolic honeycomb sandwich panels under 
uniaxial compressive loading. The test specimens were 10.00" x 
8.50", eight ply quasi-isotropic ([0,+45,90,core].) panels. The 
reinforcement consisted of either one or two additional G/Ep plies 
co-cured to the outside of eacn facesheet. Three general 
reinforcement configurations were considered: round, square and 
strips parallel to the applied load. The analytical results 
demonstrated that small amounts of reinforcement could greatly 
increase the strength-to-weight ratio. The indication was tnat 
concentrating the reinforcement close to the cutout yielded the 
greatest decrease in stress concentration. A program of 
experimental validation of the analytical results experienced some 
problems with premature panel failure caused by the facesheets 
separating from the core. It generally confirmed the analytical 
results, however. Further experimental tests on promising 
reinforcement configurations are justified based on these results. 
Properly designed reinforcement around cutouts in composite panels 
can significantly reduce the stress concentration and holds the 


promise of far lighter and stronger aerospace structures. 
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I. INTRODUCTION 


The ratio of strength to weight is one of the principal means 
of determining the efficiency of a structure. In aerospace appli- 
cations this comparison can be the most meaningful measure. The 
dilemma of designing an airframe for both strength and lightness 
has been with us since the days of daVinci. The quest for ever 
higher ratios of strength-to-weight has led to the development and 
use of high modulus advanced composite materials, principally 
graphite or carbon fibers bonded together in a polymer matrix, in 
the place of metal. 

Major airframe structural components such as wings or bulk- 
heads require cutouts for bolted or riveted attachment, access to 
interior space and passage of control and fuel lines. Timoshenko 
and Goodier [Ref. 1: pp. 78-84], among otners, point out that sucn 
holes in load-bearing structures act to greatly increase the local 
stress and to reduce ultimate strength. This characteristic is 
referred to as the stress (or strain) concentration factor (SCF or 
K). It seems the SCF may have several definitions, depending on 
the material and the researcher. In this report it shall be 
defined as the highest plane strain existing around a cutout 
divided by the far-field strain; generally called the gross SCF or 
Kgross’ Taking into consideration Saint-Venant's principle, the 
far-field strain is assumed to be equal to the strain which would 
exist in an ideal, thin, stressed infinite plate if a cutout was 
not present. Stress and strain concentrations, while inextricably 
linked in elastic materials, are not the same. However, since in 


the application discussed here, there is little appreciable 


E 


numerical difference, the term "SCF" will be used to indicate 
either the stress or strain concentration factor. 

When holes or cutouts are necessary in a structural component, 
airframe designers generally have the choice of accepting either a 
significantly lower ultimate load or greatly increasing the compo- 
nent's strength, and thus its weight. In either case the ratio of 
strength-to-weight is reduced in proportion to the highest SCF 
existing within the member.  Properly designed ductile metal 
structures mitigate the effects of SCF by plastically deforming 
under high load conditions. This response delays ultimate 
failure, but can also lead to unacceptable reductions in both 
stiffness and fatigue life. 

The metals used in aircraft construction, principally aluminum 
and titanium, can almost always be considered isotropic (many 
manufacturing processes, however, introduce some minor directional 
properties). The magnitude of the orthogonal strains (X, Y and 
shear) existing at a point ina plane isotropic panel result from 
the orthogonal stress resultant at that point and are in propor- 
tion determined by the elastic modulus and Poisson's ratio. An 
applied in-plane stress on an isotropic plate will not induce 
curvature other than, of course, the possibility of the plate 
buckling under compression. Composite plates are termed "quasi- 
isotropic" when they are composed of anisotropic or orthoceras 
laminae stacked with the directional properties arranged in a 
manner to react identically to a true isotropic material to both 
moments and inplane loads. 

Composite laminates typically lack the ductility of metals. 
The high-modulus graphite/epoxy (G/Ep) fibers in general use in 


the aerospace industry allow approximately 1$ strain (10,0004 ) 
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in tension and compression to complete failure. Depending on the 
fiber orientation, panels constructed of laminated advanced 
composites with notches or cutouts can demonstrate from sligntly 
less to much more sensitivity to holes or cutouts than otherwise 
identical isotropic metal panels. As shown by Rybicki and Hopper 
[Ref. 2: pp. 15-27], among others, this sensitivity principally 
depends on the type of weave and the orientation of the plies in 
the laminate; that is, it depends on the degree of orthotropy. 
The inherently brittle nature of advanced composite materials, 
their characteristically low strain to failure, coupled with 
manufacturing limitations make their design a far more demanding 
task than that for metals. Other characteristics, however, 
including fatigue and corrosion resistance, light weight, and 
easily tailored directional properties make the design of 
composite structures very attractive, particularly to the 


aerospace designer. 


A. OBJECTIVES AND SCOPE 

This study was designed to investigate the effect of relative- 
ly simple co-cured reinforcement of a cutout on the strain field 
and failure behavior in G/Ep honeycomb sandwich panels subjected 
to uniaxial compression. Honeycomb construction allows very light 
yet exceptionally stiff structures. The objective was to 
determine if a simple and inexpensive reinforcement geometry using 
small volumes of co-cured G/Ep lamina near the cutout could 
Significantly reduce local stress concentrations and increase the 
ultimate failure strength in the honeycomb laminate. The idea of 
local reinforcement around holes is not new; Timoshenko noted 
[Вег. 1: р. 82] that "reinforcing rings" could decrease the SCF in 


plates with cutouts. The point was to examine the reaction of an 
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advanced composite, a material whose characteristics differ 
markedly from those Timoshenko addressed. 

The research was undertaken with the manufacturer principally 
in mind. Complex or exceptionally thick reinforcement geometries 
are difficult and expensive to manufacture cost-effectively or 
with a high degree of quality assurance. This research used only 
very thin (maximum thickness: 0.028") ply reinforcement in three 
relatively simple geometries. Since facesheets with reinforcing 
olies on both sides would require machining a precise shallow 
depression in both the face of the the honeycomb core and the 
surface of the layup plate, each difficult and expensive tasxs, 
reinforcement was restricted to the outside surface of each 
facesheet. 

This study was limited to one panel size (10.00" x 8.50"), a 
single loading condition (uniaxial compression) and three rela- 
tively simple reinforcement geometries. The 1.00 incn diameter 
circular cutout was reinforced with concentric co-cured round and 
square G/Ep plies around the hole and stiffening strips displaced 
0.50 inch (1 hole radius) laterally from the cutout edge. The 
total amount of reinforcement used varied from 1 to 5 times the 
G/Eo removed from the cutout. Reinforcement was either one or two 
plies symmetrically applied to the outside of both facesheets of 
the panel. A honeycomb core was used, as it would be in an actual 
application, to increase the panel bending stiffness and thus 
eliminate the buckling of the whole panel as a mode of failure. 

The basic panel facesheets were four layer [0,+45,90] 
HMF330C/34 G/Ep cured to a thickness of 0.056 inch. Cured sheets 
were bonded to both sides of a 0.50 inch thick fiberglass/phenolic 


honeycomb core using 3M, Inc.'s AF-126 (250?F) cured adhesive. 
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The result was a very light, thin quasi-isotropic laminate: 
[07*45,90,core]. with a great resistance to bending. The 
HMF330C/34 is a woven, high-temperature epoxy (350?F) G/Ep fabric 
manufactured by the Fiberite Corporation of Winona, MN. In order 
to reduce the number of design variables the principal axis of tne 
reinforcement plies was oriented only in the direction of the 
applied compressive load. This theoretically gives tne highest 


Stress concentration and could be considered the worst case. 


B. REVIEW OF LITERATURE 
1. Background and Historical Research 

The subject of notch-induced stress concentrations in 
plates has been extensively documented. The effects reinforcement 
have on the SCF in plates have received considerably less atten- 
tion. Early research concentrated on metals (isotropic materials) 
and focused on defining the stress and strain fields around 
Circular and elliptic cutouts. Recent research has been primarily 
in characterizing the response of orthotropic and anisotropic 
materials. 

Kirsch [Ref. 3] is commonly cited as the first to 
determine exactly the stress concentration factor of a cutout in 
an isotropic material from the theory of elasticity. Howland 
[Ref. 4] applied the solution to Airy's equation in polar 
coordinates to determine the magnitude af the SCF. One of the 
earliest papers addressing reinforced holes was by Levy, Woolley, 
and Kroll [Ref. 5]. They investigated the effect of both 
reinforced and unreinforced holes on the buckling strength of 
square isotropic plates. They determined that presence of a hole 
caused only a relatively minor reduction in the buckling 


Gul C Imate) load. 
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A thorough theoretical, closed-form mathematical treatment 
of anisotropic materials with stress concentrations can be found 
in tne work of two Russian applied mathematicians, S.G. Lekhnit- 
skii and G.N. Savin.  Lekhnitskii [Ref. 6] principally addressed 
the distribution of stress around the edge of variously shaped 
cutouts in unreinforced anisotropic plates and shells under a 
variety of loading conditions. He determined that a plate with 
high anisotropy, as found in strictly unidirectional fiber con- 
struction, could produce a stress concentration factor near 9 when 
the load was parallel (0?) to the principal fiber direction, and 
slightly more than 2 with the load perpendicular (90°) to it. It 
must be pointed out that, in composite materials, the SCF may not 
have a exactly proportional effect on the reduction in the ulti- 
mate strength of the plate. Due to the composite's ability to 
redirect the load path once fibers are broken or lose stiffness 
through matrix degradation, the ultiimate strength 1s not degraded 
as much as would be expected by the presence of the stress concen- 
tration. This phenomenon is discussed in more detail in section 
11.0.4. 

Savin [Ref. 7] treated the stress and strain fields imga 
plate resulting from a cutout. He addressed the SCF as a function 
of a plate's linear material properties, ply orientation and 
Stacking, and its loading. Hole size, reinforcement and geometry 
were not addressed. A computer program was developed by Garbo and 
Ogonowski of McDonnell-Douglas LRef. 8: Vol. 3] which computes the 
stress and strain field around a cutout based on LeKhnitskii's and 
Savin's analyses. It was modified by the author for the IBM 370 


and is listed IN Appendices. 
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Substantial research in stress concentrations ln composite 
plates was done by Greszczuk [Ref. 9]. He developed a theoretical 
solution for failure stress and stress concentrations in botn 
orthotropic and anisotropic material under tension. His method 
was based on the Hencky-Von Mises distortion energy method, and 
gave both magnitude and locations of the ultimate stress.  Rybicki 
and Schmueser [Ref. 10] investigated the effects of laminate 
stacking sequence, lay-up angles, fabrication temperature and 
thickness on panel stress concentrations using finite element 
analysis. 

There is relatively little research into the effects of 
reinforcement around holes in composite plates which has been 
reported in the open literature. Virtually nothing is available 
on the behavior of notched reinforced plates in compression or on 
the effect on the type of failure of using honeycomb in such 
ک5‎ ط٥0‎ 010765 . 

Kocher and Cross [Ref. ll] demonstrated experimentally 
that titanium, graphite and steel reinforcement around a circular 
“cutout in a composite plate could reduce the SCF and increase the 
ultimate failure load in tension. Their results, however, were 
based on relatively complex, thick reinforcement geometries tnat 
have not found acceptance in aeronautical design. 

A novel cutout reinforcement metnod using bonded hoop- 
wound G/Ep disks was addressed by McKenzie [Ref. 12]. The disks 
were used to reinforce both aluminum and G/Ep plates under tensile 
loads. The method proved effective both in reducing the stress 
around the cutout and increasing the plates' strength. 

The team of Daniel, Rowlands and Whiteside [Refs. 13-17] 


did extensive experimental work in characterizing the effects of 


2] 


cutouts on a varlety o£ composite materials. They did some 
limited testing of reinforced specimens in tension and found 
proportionately reduced SCFs and increased strengtn. They 
determined that interlaminar deformations occurred in the boundary 
region of the cutout, an area they defined as extending about one 
laminate thickness from the free edge. This deformation was 

very nonlinear and could cause delamination at relatively low 


loads. Strain levels next to the cutout, prior to failure, were 





found to be higher than the ultimate failure strain of panels | 
wlthout cutouts. Based on that, they determined that the SCF did 
not necessarily produce a proportional reduction in strength. 
They recommended keeping the reinforcement close to the hole, 
using stepped diameter plies ("wedding cake") to facilitate the 
load transfer, and using 45° plies in the reinforcement where 
possible. 

Knauss, Starnes, Henneke [ Ref. 18] tested unreinforced 
0.15 and 0.24 inch thick T300/5208 panels in compression for 
unbuckled and postbuckled strength. They found that under high, 
but less than normally ultimate stress levels, the laminate around 
the hole could buckle locally, delaminate and initiate total panel 
failure. A micro-mechanical failure mode was postulated where 
limited fiver buckling at the point of stress concentration ٦ 
py local imperfections such as voids, matrix cracking or poor 
fiber-matrix bonds led to total failure. 

2. Summary Of Recent Related Rescue 

This report is the fourth in a series of investigations on 
the character of stress concentrations in composite plates made of 
HMF330C/34 G/Ep. This particular material was chosen because of 


its current use in both the Trident submarine launched ballistic 
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missile (SLBM) and the prototype Lear Fan propjet aircraft and the 
fact that it has a relatively small data base compared to other 
G/Ep prepreg material currently in aerospace use. This research 
was funded by the Strategic Systems Project Office of the Naval 
Sea Systems Command and greatly assisted by Lockheed Missiles and 
Space Co. (LMSC), Sunnyvale, CA. 

The initial project was undertaken by Herman [Ref. 19] who 
investigated the pre- and postbuckled strength of HMF330C/34 
panels loaded strictly in shear. He used a molded-in 45° flange 
around the cutout to add strength to the shear web. He determined 
that this reinforcement method was well-suited to adding stiffness 
to panels that were not buckled but that the panels did not see a 
Significant increase in ultimate strength once buckling had 
occurred. 

O'Neill [Ref. 20] demonstrated that reinforcement of only 
one face of a notched panel under tensile loading provided limited 
additional strength. Initially, the reinforcement of only one 
side of a cutout was considered attractive since only a small 
additional manufacturing effort was required. Asymmetric rein- 
forcement, however, displaces the midplane of the laminate (under 
the reinforced area) toward the reinforced side. Uniaxial tension 
tends to pull this local midplane toward the load line, causing 
out-of-plane bending at the hole, which results in high shear 
stress between plies, delamination and premature failure. The 
delamination counteracts most of the decrease in stress concentra- 
tion provided by the reinforcement. 

Paekett and Sullivan [Ref. 21] and Bank, et al. [Ref. 22] 
continued O'Neill's research examining tension panels with 


symmetric reinforcement. They showed that suitably designed 
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reinforcement which was symmetric along the axis extending through 
the panel's thickness could both reduce strain concentrations and 
proportionately increase the ultimate strength. No delamination 
was noted in their test panels. 

The work reported here extends the idea of symmetric hole 


reinforcement to compression specimens with a honeycomd core. 
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II. APPROACH TO THE PROBLEM 


A thorough investigation into the effect of reinforcement 
around stress concentrations in composite plates must examine 
various materials, hole sizes, panel and reinforcement layups and 
geometries as well as the means and directions of load applica- 
tion. This research addressed only a small portion of the total 
problem. The material, hole and panel size, layup and loading 
method remained constant--only the amount and the shape of the 
reinforcement was varied. Reducing the number of design variables 
to only two allowed an analysis of the sensitivity of the SCF to 


certain thin reinforcement geometries. 


A. METHOD OF INVESTIGATION 

To investigate the effects of co-cured reinforcement around 
cutouts, linear elastic finite element analysis (LEFEA) was 
employed to determine the strain field in each panel configura- 
tion. Plots were drawn of the strains existing on a line from the 
point of highest stress concentration at the cutout across the 
middle of the panel (the X axis) and around the cutout and con- 
tours of the three inplane strain fields (Y, X, and shear). These 
are included as figures in the appendices for each geometry. 
Specimens of each configuration were manufactured, instrumented 
with strain gages and finally loaded in compression to failure. 
The analytical and experimental results were compared and the 
failure mode of each panel evaluated. 

Total failure, in this report, is assumed to be facesheet 


delamination, separation and buckling with massive fiber failure 
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such that the panel could not withstand a full reversal of the 
load. Partial failure was facesheet delamination and separation 
witnout tne massive fiber failure and infers that there could be 


Significant tensile strength remaining. 


B. COORDINATE SYSTEM 

There are several different right-hand coordinate systems that 
have been used in the analysis of laminated materials. Analysis 
is, of course, independent of the the system used, but more than 
one student has lost his way attempting to compare methods or 
results expressed in different systems by various recognized 
authorities in the field. 

The data presented in this report is based in a cartesian 
system with the plane of the laminated panel aligned in the X-Y 
plane. Individual ply orientations are considered to be rotated 
counter-clockwise from the X axis an angle of theta (O) дедгеез. 
These plies in the layup are assigned a local orthogonal coordi- 
nate system designated the 1-2 axes. The 1 axis, also referred to 
as the principal axis, is considered to be in the fiber direction 
with the highest elastic modulus (Ej). 

Figure 2.1 shows the upper right quadrant of a typical panel 
in the X-Y (global) coordinate system as well as the ply 1-2 
(local) coordinates. This coordinate system was used by R.M. 
Jones in Mechanics of Composite Materials and Ashton Нашр ее 
Petit [Ref. 23}. Tsai and Hahn [Ref. 24] chose instead to fix tne 
X-Y axes to the ply and the 1-2 axes to the panel. The principal 
researcners in the field do not use the same system. 

The panel is oriented so that the area of greatest interest, a 
horizontal plane bisecting the circular cutout, 15 aelieneds aes 


tne X axis, where y = 0.0". The origin is assigned to the center 
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of the circular cutout. The compressive load is applied to the 
panel 90° to this plane, parallel to the Y axis, and referred to 
as Qn: The Z axis is centered at the midplane of the pane ٠٦ 
extends through the thickness toward the viewer, completing a 


right-hand coordinate system. 


C. SELECTION OF TEST SPECIMEN CONFIGURATION 

The dimensions of the test specimens were chosen to approxi- 
mate, at least in order of magnitude, a typical honeycomb panel 
with a cutout found in many aerospace applications. The overall 
size was limited by the size of the test machine and compression 
test frame. 

Hong and Crews [Ref. 25], among others, demonstrated that the 
stress concentration in orthotropic composites under uniaxial 
loading was dependent on the ratio of hole diameter to panel width 
(d/w). Whitney and Nuismer [Refs. 26 & 27] pointed out that the 
absolute hole size had a significant effect on the stress gradient 
and ultimate strength when the hole diameter (d) was less than 1.0 
inch. 

The cutout's 1.00 inch diameter was chosen, therefore, to 
limit, as much as possible, hole-size effects. The panel was then 
designed as large as practical to reduce the effect of finite 
panel dimensions and still fit into the test frame and machine. 
Hole-size and finite-width effects are addressed in more detail 
in Section II.D.3. The specimen size, 10.00" x 8.50", gavea 
diameter-to-width ratio (d/w) of 0.118 and a diameter-to-length 
ratio (d/l) of 0.100. A comparison is made in Section ITII.C.3 
between the solutions for finite and infinite plates of otherwise 


equal thickness and material constants. 
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A fabric G/Ep prepreg material was chosen because it has been 
somewhat less studied than tape and because it is finding 
increased use in airframe construction. The cured fabric laminate 
has slightly less in-plane stiffness and strength per unit thick- 
ness than uniaxial tape made from identical fibers. This is due 
еһе nature of the weave, where the fibers (or tows) are cured 
with "crimps" rather than straight. Fabric has, however, demon- 
strated significant advantages over tape in its damage tolerance 
[Ref. 28] and ease of manufacture [Ref. 29]. 

Graphite/epoxy unidirectional tape can be most effectively 
applied in flat or slightly curved structures such as wings and 
access panels. Fabric, on the other hand, lends itself to аро11- 
cations requiring high curvature or complex shapes. Tape cannot 
De used in small inside or outside radius applications without 
fiber separation, inducing matrix-rich/fiber-poor areas and 
suffering severe loss of strength. 

HMF330C/34 fabric G/Ep manufactured by Fiberite Inc. was 
chosen because it is a high modulus fabric, using Thornel T300 
graphite fibers, found in many aerospace applications. It is an 
eight harness satin (8HS) weave cloth which minimizes the number 
of fiber crimps while maintaining many of the desirable character- 
istics of cloth. Figure 2.2 illustrates some details of its weave. 

1. Panel Reinforcement Configuration 

Reinforcement of the panel cutout was of three general 
types: round, square, and strip. The round and square were 
concentric with the hole, the "stiffening" strips were centered 
0.750 incn away from the hole edge, oarallel to the applied load. 
Table I lists the panel designations, reinforcement geometries and 


amounts; Figure 2.3 shows representative configurations. 
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Figure 2.2 8 Harness Satin (8HS) Weave Cloth. 





The basic panel was a quasi-isotropic eight ply (nine 
separate layers including the core) G/Ep panel. For more simple 
comparison, the amount of reinforcement was normalized by the 
amount of G/Ep removed from the 1.00 inch diameter cutout in the 
facesheets of the unreinforced panel. The relative volume of tne 
reinforcement ply(s) was determined from this volume (0.088 in?) 
of G/Ep. Тһе round and strip reinforced panels had 5 increments 
of 100% of the removed reinforcement volume and the square rein- 
forcement had increments of 100, 300 and 500%. Тһе 200% апа 400% 


reinforcements were each two plies thick. 





Боомр 





Figure 2.3 Panel Reinforcement Configurations. 


The panel designation was devised to be somewhat descrip- 
tive of the test specimen. The first letter, P or R, refers to 
either a plain (unreinforced) or reinforced configuration, respec- 
tively. The second letter indicates the type of reinforcement: 
meme (O), round (R), square (S) or strip (H); X indicates no hole 


was present. The first numeral represents the normalized percent 
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of reinforcement, 1 to 5 for 100% to 5002 (0 india Te BR Pp 
forcement). The second numeral is the number of reinforcing plies 
on each facesheet. For example, RH42 is a reinforced panel with 
four times the removed hole volume (0.352 in?) arranged in a strip 


configuration, 2 plies thick on each facesheet. 


TABLE I 
TEST SPECIMEN MATRIX 


Panel Normalized  Ply(s) per 
Designation Type Volume ($)  Facesheet 


РХ Мо cutout or reinforcement 
2000 Мопе 0 0 
RRIT Round 100 ү 
RR22 Round 200 2 
۶1 Round 300 1 
RR42 Round 400 2 
Row Round 500 il 
К511 Square 100 1 
RS31 Square 300 1 
RSS Square 500 l 
EHI و2‎ ۲۳ 100 1 
ЕН22 SETIP 200 2 
БИ Strip 300 1 
RH42 SEED 400 2 
ВЕНЫ Бр 500 1 


— — — — = — — — — — — шш шшш шшш шш ш» шш» mm ee шш»  — — с س‎ сие чие сине — 
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Figure 2.4 shows a typical laminate cross-section from the 
midplane. Each panel was symmetric about all three axes. 
Exceptional care was required and taken during the manufacturing 
process to ensure that the reinforcement plies were placed 


directly opposite each other on the opposing facesneets. When 


38 





measured, no reinforcement was more than O.05" off center; the 


average was less than 0.02". 


D. SOME CHARACTERISTICS OF NOTCHED GRAPHITE/EPOXY PLATES 

The characteristics of composite materials differ radically 
from those of the metals they replace. As previously noted, 
composite fibers, particularly G/Ep, are by nature very brittle. 
Tensile failures in composite plates with cutouts are, almost 


without exception, load dependent. [Refs. 20-22] 


REINFORCEMENT (1 vly/facesheet) 0° 


+45° 





Figure 2.4 Typical Laminate Cross-Section. 


Compressive failure, the type dealt with here, is more 
dependent on the type and thickness of the laminate, the use of 
noneycomb to overcome the tendency to buckle, the size of cutouts 
and the presence of imperfections. The compressive failure modes 
tend to be complex, composed of one or more types of failure: 
EbLITtvy, ply delamuination, matrix cracking, etc. Stability 
failure is principally the buckling of either fibers within the 
matrix (micro-mechanical) or the structure itself (macro- 


mechanical). These test specimens and the frame were designed to 
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preclude macro-mechanical buckling (in the Euler column mode) 
since little would be learned about the reinforcement effects and 
this type of failure has been well documented beginning with Levy, 
Woolley and Kroll [Ref. 5]. 
1. Stress Concentration апе Еос s < 

It is well known that notches and cutouts in plates act as 
stress risers. For circular holes in plane elastic isotropic 
infinite plates under uniaxial tension or compression, the stress 
at the hole edge 90? to the applied load will be exactly three 
times the far-field stress. The distribution of stress around the 
hole edge and tne stress field around it can be predicted using 
Airy's stress function. Dally and Riley [Ref. 30: pp. 67-83] give 
a clear and concise derivation of the stress field equations which 
will not be repeated here. 

2. Orthotropic Effects On Stress Оа 

When an orthotropic plate with a stress riser is loaded, 
the SCF depends on the degree of orthotropy, that is, now much the 
elastic modulii change with radial direction. This is sometimes 
referred to, not always correctly, as the гана ос в оси 
subscripts "l" and "t" refer to the effective lateral and trans 
verse modulii where the lateral direction is parallel to the 
applied load and transverse is 90* to it. In the coordinate 
system used in tnis report, the load is applied parallel to the Y 
axis and the ratio is expressed as: Е ЈЕ Note that a ratio of 
1.0 does not ensure isotropy; it must be accompanied by the appro- 
priate shear modulus (Су) and Poisson's ratio (0). Гог а 
circular hole in an infinite-width plane orthotropic plate, the 
stress concentration Ky, on the cutout edge 90° to the applied 


load was given by Nuismer and Whitney [{Ref. 27: Eqn. 3] as: 
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+ 2( |E1/E, - PX) + Е1/5]. · (ОЭ) 


In an idealized infinite laminated plate, this equation 
must be equally valid in both tension and compression. This 
stress concentration factor (K ) may be considered a far better 
indication of the orthotropy of a material than the ratio E,/E,. 

The distribution of stress in the Y direction along the X 
axis (О у(х,0)) оша ожа сере frar тела) поста! scress ист) 


may be approximated using the following equation: 


0,„(х,0) = 00 1/2) (2+62+304-(к -3)(5p?-7p9)], (2.2a) 
where: 


р = а/ (х-а/2) and х > а/2. (2.2Ъ) 


The variable "d" is the diameter of the circular cutout and "x" is 
a location along the X axis (y = 0.0") when the coordinate system 
ls concentric with the hole. This relationship is a quite 
accurate polynomial approximation developed by Konish and Whitney 
Meer. 31]. 

ВЕСЕ Ot Pintte e late Width and Hole Size on SCF 

Compared to infinite plate width under uniaxial stress, 

finite plate width acts to increase the SCF. This fact becomes 
obvious in plates with a high d/w ratio. The applied stress must 
be carried by a greatly reduced net cross-section. The increase 
in SCF is due more to the net section effect than the presence of 
the cutout. Peterson [Ref. 32: pp. 110-111] gives the following 
equation to approximate the SCF at the edge of an unreinforced 


Circular cutout in a finite-width isotropic plate: 


Kp = [2 + (1-(a/w)?)1/L1-(d/w) J (2.3a) 
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This can be extended to an orthotropic plate иреге осие 


equal exactly three using: 
Ке = (К_/3) 2+01-(а4/м)31] 17+7٤0 (2.3b) 


The test specimens used in this report had a d/w ratio o£ 
О.Т. Ke was then calculated to be 3.045 for tne unreinforced, 
quasi-isotropic panel PO. At @,, = -10.0 ksi this would 
theoretically make the maximum stress -30,450 psi at the 0-2 0° 
position on the cutout (90? to the applied load) compared with 
30,000 psi predicted for an infinitely wide plate. This is an 
increase of 1.5%. Thus panel width has little more tnan a 
negligible effect on the SCF of the test specimens in this report. 

Further data that relate a plate's dimensions to its SCF 
are given by Hong and Crews [Ref. 25: pp. 8-10]. They calculated 
stress concentration factors in finite-width orthotropic plates 
under uniaxial loads using finite element analysis. They used a 
different definition of SCF, one based on the net cross-sectional 


area stress concentration 0+۵ This report uses tne SCF based 


on far-field stress or the gross SCF (Koross): The two are 
related by the equation: 
Kgross = Ket/ 1 (av) (2.4) 


То make valid comparisons with SCF data presented in this report 
selected results of Hong and Crews' analysis, converted from K,.; 


to K are listed in Table II. 


gross’ 
Their results show that quasi-isotropic layups 

(LO,+445,90]_) give results very close to the theoreticall i 7٦ 

values. Greater orthotropy in the load direction results ina 


correspondingly greater SCF. It is interesting to note that the 
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ratio of length to width (1l1/w) has an increasing effect on the SCF 
as the ratio d/w increases. 

Nuismer and Whitney [Ref. 27: p. 118] point out the effect 
of absolute hole size on panel failure: ". . . attention was 
called to a phenomenon that since became known as the 'hole size 
effect,' that is, for tension specimens containing various sized 
circular cutouts, larger holes cause greater strength reductions 
than do smaller holes." They state that the classical stress 
concentration approach does not explain such an effect and they go 
on to propose that while the stress concentration factor is the 
same, the distribution and gradient near the hole is different. 
Figure 2.5 reproduced from Ref. 27 illustrates this point. 

4. Failure Stress Criteria 

As previously noted, the SCF does not explain the "hole 
size effect" on failure. Nuismer and Whitney rejected linear 
elastic fracture mechanics (LEFM) to explain the inverse relation- 
ship between hole size and strength. They noted that while all 
circular holes in infinite width plates should have the same 
theoretical SCF, the distribution in fact changes with hole 
radius. The smaller the hole the more concentrated the stress 
near the edge appears [Ref. 27: p. 118]. Nuismer and Whitney 
proposed that when the notched stress (Q) reached an average 
value of Ou uit: the unnotched ultimate stress over some 
characteristic distance ae, that the panel's ultimate strength had 
been reached and failure resulted. This characteristic distance 
ae must be arrived at by testing a statistically significant 


number of panels. This distance ae is defined: 


ао 


1 
- О у(х,0)4х = Оу ме (2 5) 
а о О 
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TABLE II 


STRESS CONCENTRATION FACTORS-Roross (HONG & CREWS) 
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Layup EE. K L/w 0.05 0.10 0.33 О. БОИС 
[0,+45,901. 1.00 3.00 i 3.00 2.74 3.33 402 4 
[0,+45,90]. 1.00 3.00 2 3.01 3.03 3.49 ی٦‎ 
[90] 0.07 1 2.48 2.51 2.97 3.75 а 
[90] 0 4-9 ٦ 2.48 2.51 2.97 3.78 42 
pase es 100 2.06 1 2.88 2.93 3.38 3.84 m 
ЗН e 1.00 2.06 2 2.88 2.92 3.36 3.80 5.28 
8 1 ۶ 1.00 DIS 1 4.78 4.69 5.61 5.08 6.93 
о 1.00 3.78 2 4.82 4.84 5.22 >٦ 
[07 1349 Ессе 6.36 6.07 5.24 5.82 m 
[0] 7 3 s +٦ 6.44 6.44 6.54 7.3) 9: 
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х Тпалсакез ап Е./Е,. ratio of 2.00, sout now 
quasi-isotropic laminate. 





(Reproduced from Reference 27) 


Figure 2.5 Hole Size Effect on Normal Stress Distribution 
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The ratio of notched to unnotched ultimate Strength (C/O ui: 


for infinite plates is: 
gc. /0, » 2ü- E)/E2- E?- E^ +(к ,-3)( £9- E9)] (2.6a) 
where: 
š - d/(d/2 +a) and x > аә. (2.65) 


Nuismer and Labor [Ref. 33: p. 55] determined that for 
AS/3501-5 G/Ep (tape) in compression this characteristic length 
was 6.2 mm (0.24"). They also note that the characteristic length 
for tape in tension was only 2.3 mm (0.091"). Test data provided 
by LMSC indicates that for HMF330C/34 fabric G/Ep this character- 
istic length is close to 7.3 mm (0.33") in compression. 

ЕВЕ aon and interlaminar Stresses оп Failure 

Isotropic materials may be modeled using classical plate 
theory neglecting out-of-plane stresses (+z in this coordinate sys- 
tem). Orthotropic materials, however, develop» complex inter- 
laminar stress fields near the edge of a cutout. The subject has 
received much theoretical attention [Refs. 34 through 36]. Tang 
[Ref. 34: p. 1631] states that ". . . radial and shear stresses of 
each layer along the contour of the hole are in general not zero 
because there exists a three-dimensional state of stress at the 
free edge of each layer which the plane stress solution cannot 
predict."  Greszczuk [Ref. 9: p. 372] pointed out that "In ortho- 
tropic and anisotropic plates containing openings, the failure 
will take place not as a result of stress concentration, but 
rather as a result of interaction of various stress components." 

Under uniaxial compressive loading the laminate will have 
a Poisson expansion induced out-of-plane tensile stress со) 


which is highest at the hole's edge at point of the greatest 
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stress concentration. This stress is added to anv local stress 
due to machining and imperfections and combined tend to hasten 
delamination and the ultimate failure. In the experimental 
results reported here it was not possible to effectively quantify 


the effect on failure of this out-of-plane stress. 
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III. COMPUTATIONAL ANALYSIS 


Before an experimental program could be developed, it was 
necessary to understand and be able to analyze the strain field 
resulting from a cutout in a representative panel and to be able 
to predict the reaction of test specimens to an applied com- 
pressive load. Three analysis methods were used: classical 
laminate theory, tne linear elastic stress function and linear 
elastic finite element analysis (LEFEA). 

Laminate analysis provides the basic stress-strain relation at 
a point, once the material properties of each constituent ply are 
Specified. The stress function was used to predict the theoreti- 
cal stress-strain fields in an infinite unreinforced orthotropic 
elastic plate with a circular cutout. These two can be solved in 
closed-form and require relatively little computation time using 
modern computers. The finite element method allows detailed 
analysis of reinforced finite-width reinforced panels, but 
requires a significant allocation of computer resources for an 
accurate representation of the strain field. 

Tnere are several coordinate systems and notations in general 
use in laminate analysis. The following section presents the 
method used in this report, explicitly defines the notation and 


Gives justification for some of the assumptions that were made. 


А. LAMINATE THEORY AND ANALYSIS 
Laminate theory seeks to predict the properties of a multi- 
directional composite laminate based on the properties and 


E-Hntation of its constituent lamina. Individual laminae are 
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usually either unidirectional (tape) or woven (cloth) fibers 
embedded in a polymer matrix (generally a thermoset resin whose 
molecules are linked 1п three dimensions and which exhibits 
elastic properties in normal use) and tend to have strongly 
directional properties. The theory assumes that the state of 
stress is plane, displacements are small compared to laminate 
thickness and that strain is much smaller than unity. 

Pipes [Ref. 37: pp. 4-1, 5-1] presented the micro- and macro- 
mechanical models that are the basis of the theory. An anisotro- 
ріс material's elastic response at a point to applied stresses may 
be defined using generalized Hooke's Law. The constitutive rela- 
tion 1S Equation 3.1, where e and € L1 are the components of 
the stress and strain tensors and Сізкі is the tensor of elastic 
constants. Using this most general of equations there are 37 or 


8l material constants. 


3 
ва > 2, Cisi Ека (3.1) 


3 
Кі 1-1 


This equation may be greatly simplified using the symmetry of 
stress and strain and the requirement that the strain energy 
density function be positive definite [Love, Ref. 38: pp. 97-111 
and Feynman, Ref. 39: v. 2, ch. 31-7] reducing the independent 
elastic constants from 8l to 21. Symmetry reduces both GE and 
€; from nine to six different values. Feynman explains Ca sS 
elastic response of a crystal with no symmetry an e ы 
can be completely defined using 21 independent coefficients. The 
notation can be contracted using the following convention where 


the inde: ш = 6 
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= О. O23 = 04 о 05 012 = 0 


ЕЕ є; 6.3 - €, € 14 = Е € . € 


(352) 


The constitutive relation can ncw be expressed as: 


6 
Ug - >. сіз 6. (3.3) 
j=l 
This is a sixth-order symmetric matrix (where Сі; = Cji). 
Idealized thin laminate theory neglects stress and strain in 
the +z direction the equations are reduced to plane strain and 
stress, further contracting the elastic constant tensor toa 
third-order symmetric matrix. In orthotropic systems (axes at 
right angles to each other) the "1" direction is the principal 
fiber direction (or the direction with the highest elastic modu- 


lus), "2" is 90° to it and "6" is the shear in the 1-2 plane. 
01 =. Q: . Е. е (3.4) 


The matrix [Q] is termed the reduced laminate stiffness matrix, is 


symmetric and is related to [C] by: 


С. С. 
moch _— 13 —J3_ (i,j = 1,2,6) (5.5) 
С 
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The matrix [Q] may be expressed explicitly in terms of modulii 
and Poisson's ratios: 
h El oM 


О e ro e даг" 5)‏ سح کلک کا ا ا اہ 


es >} on И, 
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m Е2 95)  ء‎ : 0 "n. 

29 = С б MEE . 
7اط‎ 21 Отона ОСБ 

Св Sip Че = “Glee 


To determine the ply's elastic response defined in the 
laminate coordinate system (X-Y), both the stress and strain 
vectors must be rotated an angle @ about the "3" axis (Note: the 


"3" and Z axes are colocated): 


/ 
О, ғ x. Taig О. (1 = х,у,ху) ЕР 
1- 2 6 
! 
6: E 3 Teij Е. al УИ ХЫ) «ос! 


ј= 1,2,6 


The transform matrix [Ty J (for the case of stress) 1s derived m 


the trigonometric relations: 


д2 n? 2mn 
2 2 m = cos 
Toij = n m – 2тап. (3.9) 
22 n= sin 
-mn mn m^-n 


Recall tnat engineering shear strain (Y D drci T EE 
tensorial shear Strain (Ес) Бу абе Е € < = 2. The 
strain transform matrix elements Т.с апа Т... ресоте тп апа Тесу 
апа Те-2 become 2mn. Using matrix algebra, the now transtormed 
reduced laminate stiffness matrix [Q'] can be expressed in matrix 


form as: 
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трол Тт (3.10) 


A laminate is built from the stacking of a number of these rotated 
plies. The designer may easily tailor the laminate using various 
ply thicknesses and orientations. 

The integration of each ply's [Q'] matrix tnrough the laminate 


thickness (h) gives the normalized inplane stiffness matrix: 


h/2 
[А] = Кол а2 (ОО) 
-h/2 


Stress and moment resultants are defined by integrating stress 


through the laminate thickness: 


32 

IN} = od (ш?) 
-h/2 
h/2 

|М| = 1O} ZZ. (3.13) 
-h/2 


The stress resultant vector INI 15 related to the strain 
vector {€} by the laminate inplane stiffness matrix [A] in 


equation 3.14: 


Ny МЕ = 
Nxy Ад ^62 Ае € xy 


The laminate, while thin, demonstrates resistance to bending 


governed by the ply stiffness and the square of the distance from 
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tne midplane (+2)2. Integrating throuçh the laminate's thickness: 


1052 
[D] = [Q'] z^ az. (З.Д 
о 


Laminates with unsymmetric layups (where opposing plies at +z 
do not have identical thickness, properties and principal axis 
orientation) exhibit coupling between strain and curvature (Xx). 
This follows since each side of the midplane exhibits different 
material properties. Any applied inplane stress will induce some 


curvature. The bending-extension coupling matrix [B] is: 


n/2 
в] = ТО dz (3.16 
-1/2 


It follows, therefore, that in perfectly symmetric laminates [B] 
must evaluate to zero. 

The combined pending-extension properties of a laminated plate 
can be expressed as a sixth-order symmetric matrix which relates 


Stress and moment resultants to strain апа спкларвцие- 


(3.179 


1. Laminate Properties 
LMSC provided the initial data on material properties of 
cured HMF330C/34 G/Ep fabric. In order to validate it for thus 
program, a solid panel (PX), one without the 1.00 inch cutout, 


was manufactured and tested. The laminate material properties 
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required slight revision (less than 42$) to match the the actual 
response of the solid panels to loading. These results are 
discussed in detail in Section V.B.l. Тһе [0,+45,90,соге в solid 
laminate exhibited different modulii in tension and compression. 
In addition, it exhibited a slightly nonlinear stress-strain curve 
in compression (see Table VI and Figure 5.4). The elastic modulus 
parallel to the applied load (principal modulus, E) varied from 
7.8 to 6.5*109 psi as the applied load varied from O to panel 
failure at -57 ksi; as the load increased the panel stiffness 
monotonically decreased. This characteristic is most probably due 
to the woven plies (Figure 2.2) compressing within the elastic 
matrix, but it was not further investigated. 

The finite elements chosen for this analysis assumed 
linear elastic material properties. Nonlinear analysis was 
possible using different elements, but would have yielded little 
more accuracy at a tremendous increase in computation time. At an 
applied far-field stress (ен 90 ksi the stress induced in 
an unreinforced quasi-isotropic panel with a cutout varies from 
EIU to -30 ksi and thus E, would vary from 7.46 to about 6.95 msi. 
Since the compressive stress field and thus the material proper- 
ties vary continuously over a panel with a cutout, it became 
necessary to select one principal modulus, indeed, all the 
material constants (Ei, Ер, G15, P2 and V) for Use in he 
FEA. The material properties listed in Table III are valid (at 
70°F) throughout the range of tension but in compression they are 
exact only at -15 ksi; for other values they are approximate but 
introduce only a small error. 

Jones [Ref. 40: pp. 16-21] discussed the bimodulus 


phenomenon and proposed an improved analysis method he called tne 
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TABLE 111 


MATERIAL PROPERTIES OF HMF330C/34 CLOTH G/EP 


Tension Ei: 10.9x109 psi E>: 10.3x109 psi 
Compression Е]: 10.2х106 psi E5: 9.6x109 psi 
Shear 612: 1.0x10° psi 
Poisson ratio Vio: UO? V: 0. 09 
Tnickness t: 0. 014 inch (fully Cured SUS 
TABLE IV 
LAMINATE STRESS RESULTANT AND MOMENT PROPERTIES 
(COMPRESSION) 
A MATRIX 
8.876E-T05 Евы 1.362Е-01 
757 8.876E U5 -1.467Е-01 
1.362Е-01 -1.467Е-01 3:3349 
B MATRIX 
6.250E-02 0.0 1.221Е-04 
O0 ©.250Е-02 1. 221E-04 
1.2212-04 1.221904 5.078Е-02 
D MATRIX 
6.90 7Е+04 l-ZLLETO4 237 
1. 711E+04 6.808E+04 3. 2068701 
>. 2 РОБОМ 3926 9E 2.580Е-04 
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weighted compliance matrix. If a more complete analysis is 
required, this model should be considered. 

The bending-extension matrices (Eqn. 3.17) were calculated 
using conventional thin laminate analysis (based on experimentally 
derived material properties) for the HMF330C/34 [0,+45,90,core в 
К але in compression. The results (the [A], [B] and [D] 
matrices) are listed in Table IV. 

The symmetry of the basic laminate is apparent from the 
magnitude of the [B] matrix particularly in relation to [A] and 
[D]. The reinforced laminate also had [B] = O since it was 
symmetric. The very small relative values of the elements of the 
LB] matrix (as well as elements A3]: A45, Ау; апа A54) are more 
due to round-off error in the computer, using single precision 
numbers, than an indication of an unsymmetric layup. 

Pipes [Ref. 37: pp. 5-4] notes that when analyzing com- 
posite laminates it is often more convenient co treat them as 
homogeneous plates. For symmetric laminates it is possible to 
express orthotropic material constants in terms of the inplane 
stiffness matrix [A]. The laminate material properties may be 


determined in the X-Y plane from [A] using equations 3.18 through 


ER. 
Е Ал (пао) (3.18) 
Е, = (А11%А3; - А, 0) (3.19) 
Vy = А11/%52 (3720) 
زس فا‎ = Ај2/51] (82221) 
s hi (3865) 
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Table V lists the (experimentally derived) panel material 
properties at -15.0 Кз1. For the purpose of linear elastic 
analysis these are assumed to be constant over tne stress field 
for the particular laminate at any load. When these modulii were 
used in the finite element analysis (Он = -10.0 ksi) the maximum 


error in strain at any point in the field was less than ۰ 


B. LINEAR ELASTIC STRESS FIELD SOLUTION 

Savin [Ref 7: Chapt. II] gives a solution for the stress 
distribution in various anisotropic plates and beams with cutouts. 
Garbo and Ogonowski [Refs. 8 and 41] coded the 00.۶ іп 
FORTRAN for the case of a for the case of a thin, infinite-width 
orthotropic plate with a circular cutout. Their progran; revo 


by the author for the IBM 370, is listed in Appendix P. 


TABLE V 


MATERIAL PROPERTIES OF THE LAMINATES 


mm — аа шы шы чш — —  — — — — Y — — س ا ا‎ Р س س‎ — ce ee ee ہت سس ہی‎ = с о о шь [l —— 
—— صا — — — — — — — — — — — — — — — — — — — صا‎ д =— — — — — — — — — — — — — — — ee — — — — — — = = — ee 


Layup Plies E. Ey Gxy V ху SC 
[0,+45,90,с]. 8 220 26 25 0.321 309 
[02,+45,90,с]. 10 7.94 LED 2.40 0.269 Зе 
LO а US ge L2 8-6 Sl 25187 0.236 3. 


(Modulii *109 psi) 


сш чт ш: шж с=ш=. = — — — — — — — ee ==>‏ ج — — — — — чш шш ЫР СР‏ — — — — — — — ہت ہس ہے ше чш‏ — ص س س ب ص س 
— — — — — — — — — — = — — — — — — — — — — — ا сә» сш — — — — шт —— — — шш — — — — и‏ > == = = >= = = ہج ш ш® ст‏ — — 


Note: The 0.50" thick honeycomb core (c) had no effect on the 
inplane modulii. The panel nad an 8 ply layup except under 
the reinforcement. The 10 and 12 ply layup gives the mater- 
ial properties under the one and two ply reinforcement. 

The general biharmonic equation for an orthotropic material is 


given in Equation 3.24. The S coefficients are members of tne 


third order laminate compliance matrix [S], the inverse of [A]: 


{ef = [s] {N} 255 


5 7ج‎ DE 


922 дх4 ےت‎ 2526 9х3 ду E o дх дү? 


ДЕ ДЕ 


= Caan T -A = 0 а) 


Based on the original research by Savin [Ref. 7], Garbo and 


= 


Ogonowski point out that the stress function F depends upon the 


roots of the associated characteristic equation: 


ШІ) апа EF5(Z5) are the analytic functions of the complex vari- 
ables 2j; =X + R Y and Z5 - X * R54Y where R] and Ry are the 
complex roots of the characteristic equation. The expressions for 


the three inplane stresses are: 


0, - 2веГв? ф1(21) 9 R$95(25)] NC) 
۱تت‎ А (3.27) 


The functions $4 (24) and ф,(2,) are defined: 


Е(25) 
Ф022) = JE (0-29) 


92. 


3 


Ф (2) + 
[T д2; 


These equations have been slightly modified from their 


original form in order to apply in this report where there is no 


internal load on the hole. For a full development and explanation 
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of the equations, the reader should refer to Garbo [Ref. 41: p. 
586 | ог зауіп | Ret. 74. 


С. FINITE ELEMENT ANALYSIS 

Finite-width and reinforcement effects cannot be addressed 
using these two methods because of the discontinuities in thick- 
ness and material properties at the edge of the reinforcement. 
Finite element analysis has demonstrated its ability to accurately 
analyze the majority oft problems in elasticity. The quality of 
the solution is, however, dependent on the size of both the 
available computer core memory and the analysis budget, since the 
quality and cost of the solution are functions of the fineness of 
tne element mesh. The solution time and cost increase at least 
with the cube of the degrees of nodal freedom (DOF) in the model 
[Ref. 42: pp. 391-402]. Іп а full three-dimensional analysis each 
element node point may be displaced in the X, Y and Z directions 
and also rotated about each of the three axes. Thus there are six 
possible DOF per node: three displacements and three rotations. 
The dimension of tne stiffness matrix is the sum of the degrees of 
freedom at each node point in the model. 

The structural finite element analysis method, in its simplest 
form, is the determination of the relationship between the load on 
and the displacements in a body. The two are related by tne 
stiffness characteristics of the body. The body is divided into a 
number of smaller volumes (or areas) termed elements, each element 
is then assigned a local "stiffness" and these are then combined 
in matrix form to establish their inter-relation. The result is 
termed the stiffness matrix [K]. Each element is made up of nodes 
at its corner points which can be fixed, at which a force сай SS 


and which can deflect if not fixed. The vector of forces IF} 


acting on each node equals the product of [K] and the vector of 


deflections 161. 
fF} = [К]16] 0) 


Since the forces are generally known and it is the displace- 


ments which are desired, the stiffness matrix must be inverted: 
15} = [К ЧЕ] (531; 


The order of the [K] matrix is determined by the sum of nodal 
degrees of freedom. The matrix inversion to [K]T} is not a 
trivial computational task in any realistic finite element model. 

The art in FEA is in defining a mesh fine enough to give 
adequate solution accuracy while suppressing as many DOF as 
possible to keep the cost of solution within reason. Zienkiewicz, 
in his excellent text on the subject [Ref. 43], covers this method 
of structural analysis in some depth. 

TIT finite Element Program 

A finite element analysis program named DIAL as ме11 аз а 
significant allocation of computer time on a Digital Equipment 
Corporation VAX 11/780 was пайе available by LMSC for this re- 
search. DIAL is a flexible, general purpose finite element code 
for the analysis of two- and three-dimensional structures. It has 
a modular architecture ln which individual subprograms are exe- 
cuted as the model is being defined, the mesh generated, the 
equation bandwidth optimized and the solution found. As each 
subprogram (called a "processor") is executed, it extracts re- 
quired data from a data base, processes it, updates the solution 
and adds to the data base. This architecture provides an 


invaluaole restart capability at the last successful process which 
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te Element Mesh 


ИР 


Typical Е 


Figure 3.1 


can Significantly conserve analysis time [Ref. 44]. The following 
DIAL processors were used in the linear elastic analysis of the 
experimental test panels: 

* MESH The geometric grid of elements to be analyzed, 
called a mesh, is generated by specifying points coincident with 
the quadralateral element corner nodes in an orthogonal I-J 
coordinate system. Certain points key in the I-J system are then 
givenlocations in the X-Y plane and the MESH processor automatically 
maps appropriately shaped elements. Figure 3.1 show a typical 
element mesh. The processor allows partial meshes to be generated 
individually and then merged to each other creating a larger 
model. The heavy lines in Figure 3.1 outline these. Building a 
complete FEA model from a series of smaller partial meshes reduces 
tne manhours required to generate the model and allows more 
complex geometry. Merging adjacent partial meshes eliminates any 
redundant nodes and degrees of freedom. The panel models for this 
analysis used from three to five partial meshes. Boundary condi- 
tions are specified and DOF suppressed within the MESH processor 
to adequately simulate the structure. 

* BAND It is not necessary to store the entire finite 
element stiffness matrix; the Betti-Maxwell reciprocal theorem 
requires that tne stiffness matrix must be symmetric. ТЕ can ре 
decomposed into a lower and upper triangular matrix, recovering 
almost half the memory or storage area originally required. 
Further improvement can be gained by reordering the node numbers 
to optimize the matrix bandwidth and storing it using the 
"skyline" method. BAND offers a number of options to do this 


including Collin's and the Gibbs-Poole-Stockmeyer algorithms. 
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* SETUP The undeformed finite element data sets are 
generated and a series of error checks are done verifying the 
element grid. 

* МАТЬ The material properties of each ply are defined 
in MATL. The processor uses classical laminate theory to compute 
the bending-extensional properties. The strength of the processor 
architecture now becomes apparent. It possible to vary the 
material properties of the model without regenerating the complete 
element mesh, optimizing the band width or generating new element 
data sets. 

* LOAD The LOAD processor generates consistent load 
vectors for any combination of pressure, traction, body forces, 
inertia loads and temperature variations. It allows the variation 
of loads without regeneration of the stiffness matrix. 

* DIAL The nodal deflection analysis and stress-strain 
computation is done within the DIAL processor. It uses the total 
Lagrangian formulation method to handle geometric nonlinearities. 
FORTRAN double precision representation (64 bit) and sequential 
improvement to convergence was used to increase the accuracy of 
the solution. This insured the best possible solution but 
increased the equation solution time by a factor of about eight. 
The effect of using double precision and convergence can be seen 
in the figures in Appendices A-N where snear strain is resolved to 


as low as +0.003% of the value of E,, long the X axis. 


У 
* СЕОМ The data generated by even a small model is 

extensive and difficult to evaluate in tabular form. DIAL pro- 

vides an extensive array of post-processing choices to present the 


data in graphical or tabular form. GEOM generated the strain 


contours for each panel which are presented in the appendices as 
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well as the extrapolation of strain data from element Gauss points 
to the nodes. 
۲٣٠٠۰٠٦٠ ٦٦ ٥٦٢٢٦٢٠٢٠٢٢٠ ۶ the Finite Element Model 

Each panel reinforcement configuration required a separate 
finite element model. A modified thick-shell elastic quadra- 
lateral element was used for the analysis. It used the laminate 
material properties developed in the MATL processor and ply 
thlckness with any offset from the Z axis (specified during the 
mesh development) to define and individual element stiffness 
matrix. The greater the Z offset--the greater the resistance to 
bending. Several of these elements may be stacked through the 
thickness by merging partial meshes. Stacking meshes results in 
the direct addition of element stiffness. Element properties are 
projected to a reference plane (z = 0.0" for this model) which 
contains the nodes points. Element bending resistance is 
determined by its stiffness and offset from this reference plane. 
vantage Of this type of element is that lt allows modeling a 
thin three-dimensional laminate using a two dimensional element, 
thereby greatly reducing tne number of individual element nodes. 
The modified thick shell element's shortcoming is that it cannot 
give stresses in any of tne stacked partial meshes and the strain 
is valid only at Zz = 0. Further, thin plate theory is used which 
gives strain and the stress resultant (М) vector. Іп these models 
the true thickness varied over the surface of the panels. Plots 
of stress resultants in this case would be, at best, misleading. 
Strain was used as the basis for comparison among the panel 
configurations. 

During the experimental phase of the research, premature 


facesheet separation from the noneycomo core became an unpredicted 


failure mode of some of the specimens (panels RR22, RR42, RR51, 
8551 & RH31). This type of stability failure could not be 
predicted employing two-dimensional analysis only. It was oostu- 
lated that the [0,+45,90] facesheet layup or the 05,45,90] (or 
[03,+45,90]) reinforcement could be generating sufficient out-of- 
plane (+z) forces to tear the facesheet away from the core. 

To answer the question, a three-dimensional analysis using 
thick shell elements representing the facesheet, combined with 
isoparametric solid elements, representing the honeycomb core, for 
a 3-D analysis. These results (for panel RR22) are given in 
Apoendix O (results showed no significant out-of-plane stress). 
While DIAL can handle a 2-D element with an aspect ratio (length/ 
width) up to 20 with little loss in accuracy, an effort was inade 
to keep this ratio below three. At high aspect ratios the inter- 
polation assumptions within each element are no longer valid. The 
meshes employed were also designed to keep interior element angles 
as close to 90° as possible, again to increase the accuracy of the 
solution. Figure 3.2 shows the elements (numbered from 1 to 30) 
and node points (numbered from 1 to 117) next to the cutout. 

These elements’ dimensions remained unchanged (except for the 
added reinforcement thickness) for each model, allowing direct 
comparison among reinforcing configurations in the region near the 
cutout. Since tne plate was symmetric about all three axes only 
the upper right quadrant of each specimen was modeled. 

In the experimental fixture, the compressive load was 
applied to each 8.50" wide specimen using clamps 8.00 inches wide. 
The outside edges of each panel had 1/4 of an inch inside a slot 
in the vertical member which could not be loaded. To simulate a 


-10 ksi far-field load in the FEA model a Constant Tine fea ім) 


o4 





Figure 3.2 Elements and Nodes Next to the Cutout. 
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of -1,120 lb/in was applied to the 0.112 inchi 1 1ة‎ = 
(neglecting the core). The effect of modeling the test fixture 
can be seen ln the quarter-panel contour plots in the appendices 
(for example Figure A.4): there is an obvious stress concentration 
in the upper right-hand corner of the panel. The effects of this 
stress concentration die out rapidly as the distance into the 
panel from the line of load application increases. 
3. Interpreting Finite Element Analysis Results 

DIAL, like most finite element programs, produces 
voluminous data files giving the stress, strain and displacement 
at each element's integration points. Meaningful comparison of 
these files among the various panel configurations would be 
tedious as well as unenlightening. The items of interest were the 
distribution of stress along the X axis and around the cutout and 
the strain fields on each panel resulting from the reinforcement. 
Graphical comparison was chosen as the best method both to present 
and to compare the panels. 

The results from tests of each of the 14 test specimens 
with a cutout is presented in individual appendices (A through N). 
Each configuration has a plot of the element mesh and a comparison 
of strain both along the X axis and around the hole under a far- 
field normal stress load of -10 ksi. In addition, strain contours 
for (gr, = -10 ksi are shown for each panel's upper right quads 
and for the region near the cutout. Experimental data are 
correlated with the finite element analysis for each panel. 

A plot of a deflected element mesh is presented in Figures 
3.3 and 3.4 to illustrate the analytical and experimental оочпаагу 
conditions imposed on the test panels. Figure 3.3 shows the 


entire upper right quadrant of the panel. Figure 3.4 snows the 


elements close to the cutout. Тһе dashed lines represent the 
outlines of the elements prior to the application of the load. 
The solid lines show the elements in the panel compressed under a 
-10 ksi load. The deflections shown are, of course, an 
exaggeration of those actually present in the panel, but they are 


accurate representation in relative scale. 
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Figure 3.3 Typical Deflected Mesh Plot. 


The boundary conditions imposed on the quarter panel are 


Clear: the X axis, representing the longitudinal bisection of tne 
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panel allows no movement in the Y direction but allows Poisson 
expansion the X direction. The panel boundary on the Y axis was 


constrained in X displacement but was allowed to move vertically. 





Fiqure 3.4 Element Deflection Next to the Cutout. 


A polnt to note ls the boundary condition at the top of 
the panel. It was necessary tO firmly clamp the upper 8.00" x Í 
1.00" inch area of the panel to assure complete and even load | 
transfer into the G/Ep facesheets. The results of the boundary 
condition can be seen in the deflected mesh plot; Poisson 
expansion was not allowed where the panel was clamped. This very 


closely modeled the experimental setup. 
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In Figure 3.3 the top edge of the panel has a slight slope 
upward as X increases from 0.0 to 4.25". This is the result of 
applying a constant stress boundary condition along the edge 
rather than constant displacement, which would more closely model 
the experimental apparatus. This tends to slightly increase the 
SCF at the hole because the panel finite element model appears 
somewhat less stiff directly above the cutout than the solid 
portion. 

A test case using constant displacement boundary condi- 
tions, which closely approximated the experimental setup, produced 
less than a 0.5% increase in the SCF. Since each panel has a 
slightly different stiffness in the Y direction, it would have 
been exceptionally difficult to impose an identical load on eacn 
for comparison. 

Hong and Crews [Ref. 25: pp. 4-0] reported significant 
differences in results between constant stress and constant 
displacement boundary conditions. In the final analysis, the 
researcher, understanding the differences and the compromises, 


Inst choose the model best suited to his work. 
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IV. EXPERIMENTAL ANALYSIS 


A program of experimental verification was developed in order 
to determine if the analytical results of the finite element 
analysis represented the actual strain field. Each reinforcement 
geometry previously described was manufactured, instrumented and 


tested. 


A. TEST SPECIMEN MANUFACTURE 

The test panels were manufactured by Lockheed Missiles and 
Space Company using methods similar to those for verification of 
Trident missile structures. The initial uncured prepreg plies of 
HMF330C/34 clotn G/Ep were laid up on a stainless steel plate in a 
4 р1у ([0,+45,90]) facesheet. А precut uncured reinforcement (one 
or two plies) was then placed in position on the top of the 
uncured layup and retained in place by a small pin. Two identical 
facesheets were made for each geometry. A standard "bagging" 
process and cure cycle for the 350°F (450°K) Fiberite 934 epoxy 
prepreg was used. This included a hold in the autoclave at 360 
+10°F (455 +5°K) for two hours. Using acid digestion techniques 
this cycle typically yields a fiber volume of 62 +2% and a void 
content less than 1%. One facesheet was joined and cured toa 
0.50 inch thick Hexel fiberglass/phenolic honeycomb sheet using 3M 
Inc.'s AF-126 (250°F curing temperature) film adhesive (known as 
"Blue Glue"). Once the first facesheet and honeycomb had been 
bonded an aluminum/epoxy potting compound was poured into the 
honeycomb cells within 1.25 inches of each end. The aluminum/ 


epoxy potting compound provided dimensional stability for the 
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panel, assisted the load transfer and prevented crushing the 
honeycomb in the panel ends when they were clamped into the 
compression test frame. The second facesheet was then joined and 
the now complete rough panel put through a third and final cure 
cycle. The panel configuration, excluding the one or two ply 0° 
reinforcement, became [0,+45,90,core],. The core's elastic 
modulii in the X and Y directions were virtually nil and did not 
contribute to the panel's inplane stiffness. 

The center of the reinforcement was marked and a starter hole 
drilled with a No. 4 carbide-tipped steel drill rotating at 
approximately 2200 rpm. The hole was enlarged in steps using 0.50 
and 0.75 inch diameter carbide-tipped drill bits. The final 1.00 
inch finished hole was cut using a carbide-tipped boring head 
rotating at 1600 rom moving in depth at 0.0015" per revolution. 
This method provided a very smooth and almost perfectly circular 
cutout. Each facesheet was drilled using stiff fiberglass sheets 
as backing to minimize the breaking of fibers on the bottom ply 
when the drill bit broke through. Fiberglass tabs (8.0" x 1.0" x 
0.25") were applied on both sides at either end to provide for 
load transfer from the test frame into the panel. The rough panel 
was then cut to the specified size (8.50" x 10.00") using a 
diamond-coated circular saw. Great care was taken to both keep 
tne cutout in the center of the panel and to insure the two edges 
to be loaded were parallel. The general dimensions of the 


specimens are shown in Figure 4.1. 


B. TEST APPARATUS AND PROCEDURES 
In compression, much more so than tension, lack of attention 
to maintaining proper boundary conditions can quickly invalidate 


experimental results. Great care was taken in the design and 
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lmensions 


Figure 4.1 Compression Panel D 
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construction of the test frame to insure that it was extremely 
stiff, that the compression surfaces were parallel and that they 
would remain so during the entire compression sequence. 
ES Test Apparatus 
a. Load Application 

A Material Test System (MTS) Series 810 hydraulic test 
machine was used produce the compressive loading. The compression 
test frame was designed to be strong enough to utilize the 100,000 
lb. maximum load of the MTS machine. It consisted of a fixed 
horizontal base and vertical side posts and a sliding horizontal 
top cross member. Both horizontal members were machined from 
7075-06 aircraft-grade aluminum. The vertical posts were turned 
to a diameter of 2.000" from diameter mild steel bar stock. The 
horizontal members were fitted with a means of clamping the test 
specimens. Each had a 0.250" thick tempered tool steel base plate 
positioned to transfer the compressive load into the test specimen 
and to prevent damage to the surface of the aluminum frame. These 
load plates were carefully adjusted during installation to ensure 
that they were parallel within a tolerance of +0.0005 inch. 

А 0.014 +0.001" slot was milled in both steel vertical 
members to accommodate the panel and to allow some vertical 
movement while preventing out-of-plane deflection. The lower 
horizontal member was held fixed relative to tne frame while the 
top one was allowed to slide vertically. Bronze bushings were 
pressed into the upper and lower frame members and then machined 
to within a +0.001" tolerance. The vertical posts' ends were 
fitted into these bushings. A special effort was made during the 


design and manufacture of the test frame to keep tolerances as 


ce 


small as possible to maintain proper and repeatable test boundary 


conditions. 
Figure 4.2 shows some details of the compression test 
frame. The following numbers indicate some of the parts and 


features of the frame and correspond to the numbers in Figure 4.2: 
Tempered tool steel compression support plate. 

Bottom horizontal frame member. 

Vertical steel post. 


Slot to hold edge of the test specimen. 


e О 


Bronze bushing. 


The test frame was allowed to "float" in the MTS 
machine. Steel bearing surfaces were fitted to the top and bottom 
which allowed the test fixture to slide parallel to the floor for 
centering. These also eliminated the possibility of transfer of 
any moments from the MTS machine to the frame. Each steel bearing 
block was made of three pieces: опе threaded to mate to the MTS 
moving piston, a circular 2.000" diameter lubricated cylindrical 
bearing and one threaded to mate with the test frame. Figure 4.3 
show the test frame positioned in the MTS for a test. The bearing 


blocks can be seen in the figure between the test frame and the 


machine. 
b. Strain Measurement ат 
A Vishay Measurements Group, Inc System 4000, shown in 
Figure 4.4, was used to record the strain gage indications. 1; 


consisted of a Hewlett-Packard 9825В microcomputer linked through a 


Measurements Group, Inc., Instrument Division Model 4200 
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controller to Model 4270 strain gage scanners. Integral software 
provided for gage identification, calibration and strain reading, 
conversion and printing. The entire experimental test station is 
shown in Figure 4.5. 

2. Instrumentation Procedures 

Each panel was instrumented with a variety of strain gages 
principally located along the X axis and oriented in the Y 
direction. The primary purpose of the reinforcement was to reduce 
the maximum strain, and thus the SCF, at the edge of the hole 90° 
to the applied load. The 0.50 inch honeycomb core was used to 
eliminate panel buckling. The panel was designed to maintain, as 
closely as possible, equal strain on opposite facesneets. Тһе 
gages were located on either side of the hole, but on only one 
facesheet. In retrospect, gages on both sides of the cutout on 
botn facesheets would have given additional insight into the 
fallure mechanisms. 

The choice of strain gages was based on the strain 
gradient near the cutout, the panel strain field and the heat 
transfer properties of the G/Ep panel. 

a. Measurement of Strain Near a Cutout 

The measurement of strain near a cutout in the 
presence of very high strain gradients is not a straight-forward 
exercise. Reference 45 points out that an electrical gage 
effectively integrates the strain field under its grid. When that 
field changes very rapidly the accuracy of the measurement can be 
strongly affected. The typical strain field studied here demon- 
strated gradients as high as 16,000 microstrain per inch within 


0.025" of the cutout ae” =10:0 ksi far-rrewm Tc 
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Figure 4.5 


As recommended by Ref. 45, a number of techniques were 
used to accurately measure the strain field along the X axis. The 
smallest possible gages were chosen for use for next to the 
cutout; a series of in-line gages close to the hole gave strain 
gradient data. Special care was taken to accurately measure tne 
position of each gage. A Rockwell Corp. electronic, digital- 
readout gurney gave the gage center location to within +0.002 
inch. This resulted, at the -10.0 ksi test point, in about 0 
microstrain or +1% maximum uncertainty in strain due to gage 
position error. 

b. Strain Gage Excitation Level 

Strain gages require some electrical excitation to 
allow measurement of the change in resistance in the gage grid 
caused by tension or compression. This results in some degree of 
resistive self-heating within the grid. This self-heating charac- 
teristic can cause significant drift in indicated strain from the 
true value. when measuring strain in most metals, there is little 
heat buildup due to their superior heat transfer characteristics. 
What little there is can usually allowed for by self-temperature- 
compensation (STC) in the gage. STC requires the matching of 
coefficients of thermal expansion (a) of the gage and the 
specimen. The heat transfer characteristic of G/Ep is) lew 
compared to metals: The temperature under a gage can rise enough 
to invalidate the indicated strain reading. 

Reference 46 recommends a maximum power density of 
0.1-0.2 watt/in* for materials with low thermal conductivity such 
as G/Ep. Power density (PD, watts/in*) is a function of gage 
active grid area (A, іп2), gage resistance (R, ohms), and gage 


excitation level (E, volts) according to the relation: 
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PD = B/(4*R*A). رای‎ 


A typical 120922 , 0.040 in? gage at 5.0 volts excita- 
mon nas a PD of 1.3 watt/in^. As noted above, gages with a small 
grid area were necessary to accurately measure high-gradient 
strain. А Measurements Group, Inc.  EA-xx-030CM-030 gage (A - 
0.0025 in?) [Ref. 47: p. 7L] which could meet the size require- 
ments has a PD in excess of 21 watts/in^. 

Clearly, high strain gradients and composite materials 
require extreme care in selecting and using strain gages. A 
combination of lower than usual excitation levels and higher gage 
resistance were used too in this research, where required, to keep 
the power density within acceptable limits. 

C. Strain Gage Application 

Gages were applied to the panels one facesheet 
on either side of the cutout along the X axis. They were applied 
in accordance with the manufacturer's recommended procedures 
[Refs. 48 through 51] using M-Bond 200 adhesive. Figure 4.6 shows 
a typical strain gage layout on a test panel. It should be noted 
that, although similar, each panel had a unique gage layout. 
Several gages were mounted at points other than along the X axis 
to verify the analytical strain field. 

B Test Procedures 
Тре test specimens were allowed to age for 180 days at 
70°F and 50% relative humidity to reach hygrothermal equilibrium. 
Immediately after the strain gages were installed the panel was 
mounted in the compression test frame and loaded to failure. 
The test consisted of initially loading the panel to -2000 
Pst to set it in the test fixture. The load was then removed and 


the gage readings reset to indicate zero strain and then 


31 














N —— — ہر یں ہے ہدمع‎ 
-j — "3 = vero رفا لہ ون‎ сем ага мы سس‎ ee ee й اس‎ та 7 Ay .اڈ‎ — а" 
, ( | ' a 4. مد ا ہا کا‎ y ] کے ے 30ج‎ та sya У | = I 
| $ 
| ۸01۰ 
i rt 0) 
4 : 1) $4 
i ) 
ay : у E š 
і 0 f 
= ! А т 
© А + 1) 
"s q м . = {з 
n] № О de ° 0 0) 
о » ei > С 
t: ust "D С, 9 ^ 
و‎ Uc (U! 0) = 5 
~ тт} r4 0) +) О 
a ge? t v % С 
2 с Q m 1 
а tf) с! i 
= ad - ord ff) (5 м 
т z т} O, 0) 
Ü c С 
П б Ф O 2 m 
пос с с ш 
2 уз لد‎ О ia 
l 42 0 CN 4) 
ео tn y 
18 0 oe e E с, 
О О n d c Е 4) 
ЕЦ 
g 9 vo U O 43 
а) т 11 б Ф orf G 
- 760 00-21 U 0) 
ЕЧ ہے‎ fa ЛО, (Ú }4 
rd nj Q | 
(7) Т 
tf) (л 0) r-d OG 
0) Fu > м Эң 
٠ - ہے‎ 0) O a 
N مت‎ x 1 0) My е 
z 1) m л 9 T 0) 
TEE MU 2 
D г. E "m А, 9 j 
n V ۲ہ‎ 0) 
с в J = (Y, Q 
E Е 
4 (1) 4 
ا‎ (f) 4) EJ E 
° --1 tf) > 
uj - G (Ú С 
U o UO u Ф 
pou v О з e 
7 2 “41 r+ GE 
о 1 
О 4) 0) t Z 
E $ a 5 q تہ ا ا وا‎ i "ñ 3 1 
д th YƏ. нь 
EC: о < | £m n nort eb وی‎ = 
a 0 О 0) о MOL Tuve b ha ab OU Sp € oT re HUP QN Ege ^c Mun aene К Ж 
\! = af U) ч 


aqu 


1 


UDS 


In ats 


2 


used 
Э 


A oe 


25 


~ 


у 


typical test required 15-20 minutes to complete. There were some 
variations in this straight line load procedure which are noted in 


the appendices for the affected panels. 
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V. DISCUSSION OF RESULTS 


A. COMPARISON OF COMPUTATIONAL RESULTS 

All computational analysis were done at a far-field applied 
stress of -10.0 ksi. These analyses assumed the material had 
linear elastic properties. This assumption was adequate to 
reasonably predict the strain field below the material yield 
point. There are some significant nonlinear yield characteristics 
of composite materials that require more sophisticated treatment 
than is given in this report. 

l. Open Versus Closed-Form Analysis for an Unreintoncesme ms 

A comparison of open and closed-form strain distribution 

around a cutout in an unreinforced panel (POUQ)is shown in Figure 
5.1. The lines represent the infinite plate width strain computed 
uSing the stress function (Equation 3.24) by the FORTRAN program 
"RBSFM" in Appendix P [Refs. 8 & 40]. The triangular points 
indicate the LEFEA strains at the node points for tne finite-width 
(8.50") plate The effect of tne finite panel width and the con- 
stant stress loading boundary condition may perhaps be more easily 
seen in Figure 5.2 where the FEA strain results are represented by 


crosses. The maximum FEA computed strain is higher increased at 


the edge of the cutout (x 0.50") compared with the closed-form 
results. At distances more than 2 hole diameters away from the 
cutout (x >_ 2.0") the FEA model gives slightly less strain. The 
differences between the two analysis in Figures 5.1 and 5.2 are 
small; it is the similarity of the two that 15 stn ae. 


Tne increased FEA strain at the hole is due to the constant 


stress loading boundary condition. A constant displacement 
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Figure 5.1 Strain Comparison Around the Cutout. 
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boundary condition would have almost entirely eliminated even this 
small difference. The reduced strain toward the panel's free edge 
ls due to not applying the load to the outter 0.25 inches o£ the 
panel's top edge. 

The point of these comparisons is to validate the finite 
element analysis method, tne type element and the configuration 
chosen. It is assumed that the computational results are as valid 
for reinforced panels. 

2. Finite Element Analysis Results 

Table VI summarizes the most important data from the 
LEFEA. The three maximum strains (Y, X and shear) are given for 
each configuration as well as the finite-width stress concentra- 
tion factors. The locations at the edge of the cutout are listed 
at the bottom of the table. These values are best used as a means 
of comparison among reinforcement geometries, not for and exact 
prediction of the micro-strain at the edge of the cutout. Recall 
the assumption made that the compressive modulus was constant for 
all strain. Two SCF's are given, one for a theoretical "infinite" 
plate and one for the 8.50" width panel used in this research. 

These stress concentration factors are theoretical only. 
They are valid solely for a totally elastic strain field. Nuismer 
and Labor [Ref. 33: p. 50], among others, point out that atem и 
strains (in the case of HMF330C/34 at strain in excess of SOD D Rp 
the fibers immediately next to the hole at the SCF begin to fail 
and transfer the load througn the matrix to adjacent fibers. 
Compressive failure usually consists of matrix cracking, fiber 
micro-buckling, ply delamination and the transfer of the load from 
the failing region next to the hole away to fibers/matrix able to 


sustain tne load. Failure in this manner is difficult tora ars 
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Figure 5.2 Strain Comparison on the X Axis. 
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using linear methods because of the rapidly changing material 


properties during the process. 


TABLE VI 


FINITE ELEMENT ANALYSIS RESULTS 
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Maximum Strain Around Strain Concen- 
the Cutout tration Faetos 
Panel  ---------------------- ----------------- 
Designation eps-y ерз-х eps-Xy м = СО w = 8.5" 
РО00 -4230 1596 -3253 3.00 За 
RR11 -3739 1420 -3211 2.95 2:75 
2822 -3363 1284 -3113 2.39 2245 
0 -3005 09 -3106 2.56 2505 
RR42 -3231 1186 -2998 2.29 2437 
RR51 -3539 1207 -3056 2007 2.00 
ROLE -3719 22 -3170 2.04 2214 
Roo -3545 1258 - 3058 ge] 200019 
RS51 -3405 1193 -2999 2:49 2 ۶۶ 5 
RALL -4261 0 -3298 3.02 а 
RH22 -4097 T22 - 3188 2.31 32502 
RH31 -3983 1545 -3097 2032 2.92 
RH42 -3727 1045 -2947 2.64 2 
7151 -3997 1475 -3094 2:93 2.93 


Figure 5.3 shows a comparison of the maximum Y direction 
strains (eps-y in Table VI) at the edge of the cutout. These 
correspond to the theoretical stress concentration factors. 
Several facts become apparent: 


_ = Panels RR22 and RR42 gave thek ٣ 7/۲ 
reduction іп SCF. They had most of their reinforcement concen: 
trated next to the cutout in 2 plies (thick) per facesheet. 


" In no case was the 500% (single ply) reinforcement 
appreciably superior to the 300% (single ply) configuration in 
reducing the SCF. Reinforcement relatively far removed from the 
hole edge added little strength to the panel. 

Ш The square reinforcement configuration provided very 
slightly more strain reduction compared witn equivalent round. 
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Figure 5 
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B The strip reinforcement resulted in about 12$ higher 
strain for the same amount (percentage) of reinforcement compared 
with the other two configurations. 

B. COMPARISON OF EXPERIMENTAL RESULTS 

1. Solid Panel 

Panel PX) was tested to provide a basis for comparison 
and an indication of the ultimate strength of a panel without a 
cutout or stress concentration. The panel was subjected to two 
loading sequences: the first up to Ff, >>" ин 
thirds the estimated ultimate load) and the second to failure at 
-57.0 ksi. Two load runs were used to determine if there was any 
residual damage from the first load. After the first run, resi- 
dual damage would be indicated by reduction in inplane stiffness 
(the effective modulus E) due to matrix degradation. Table VII 
shows the results of the test. The first run was a maximum load 
of approximately 80% of the ultimate; no significant difference is 
apparent between the two runs. It appears that this G/Ep material 
is elastic, at least up to about 80% of its ultimate compressive 
strength. 

The monotonically decreasing stress-strain curve (noted in 
Section III.A.1) is significant. It most probably results® Bos 
the decreasing ability of the crimped harness weave fibers to 
carry the compressive load as the load increases. The close 
correlation in strain (Eps-y) and modulus (Ey) between runs 1 апа 
2 would seem to eliminate matrix cracking or delamination, at 
least below -45.0 ksi, as a source of the nonlinear behavior; it's 
nonlinear and elastic. 

Figure 5.4 shows the stress-strain curves for both load 
sequences. Test sequence No. 1 is almost exactly duplicated by 


No. 2. These curves came from the average of 10 gages mounted 


90 


(ut/ui) ۱۱18-۷ 
00001 0008 0009 0007 


озпива ој '2'омрвој O 
І8і4 очу ој '|'омрвој O 
QN3931 NIVHLS 


189 00049 :9зпива фечва |8101 
159 009 "95 :92 пе 19414 ејаирпу 18214 


О00Е 





SSƏ1S ел!$5э14и105 р|э14 184 SA UIB11S-0J9IlJ 


1n01n2 ом печеа шеа :QOXd I9ued 


(58 0001) $59235 aA4issa1duio2 petrddw 


Figure 5.4 Stress-Strain Curve for Solid Panel 
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transversely (on the X axis) on the panel. “he gages were all the 
various sizes and resistance values used on the other panels. The 
standard deviation of all the values was within +4% of the average 
for each load. This could be taken as the typical limit of 
accuracy for any one gage. When including consideration for the 
position error noted above it would not be unreasonable to 
consider any experimentally measured strain to be within about +43 
of the true value. It is doubtful that more accurate measurement 


is possible without taking extraordinary measures. 


TABLE VII 
PANEL PX TEST RESULTS 


Run #1 Run #2 

Load eps- Е У eps- У 
(-ksi) (ше) (ші) 1% (es (mu) 12 
5 644 ھ2‎ 6 ОРЕ 636 7.86 0.336 
Lo 1341 7240 OLS 1342 7.45 01925 
15 2065 7.26 QUEE 2060 уде 02820 
20 2799 ت70‎ +0 а 2793 "ENS 0.319 
2:5 3543 7 6 02315 3544 Taos Owe Le 
30 4323 6.94 О 4302 6.97 0.318 
35 5092 6.87 Ото 5073 6.90 0.319 
40 5879 6.80 O55117 5891 6.79 OQ. rS 
AS 6683 6.73 0.318 6653 6:75 0:920 
50 -- -- -- 7485 6.68 ھ9‎ I 
55 -- -- -- 8328 6.60 ОЛЕР 
5S7 == ЕЕ -- 8702 6.34 0.322 
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2. Panels with Stress Concent rat ons 
The results of the experimental program are summarized in 
Table VIII. Individual panel experimental and computational 
results are discussed in the appendices. The loads in (ksi) are 


listed for the first audible ply failure (FAPF) and ultimate.  Tne 


22 


FAPF was nothing more tnan the first "pop" heard during the 
loading sequence. While this hardly seems to be a rigorous defi- 
nition, in every case the FAPF appeared to be a predictor of the 
ultimate load. Stress concentration factors (SCF) were taken from 
the finite element analysis. In the strength reduction column the 
calculated value came from equation 2.6 using ae - 0.33", the 
value determined by LMSC for HMF330C/34 cloth G/Ep. Nuismer and 
Whitney [Ref. 26: pp. 122-3] state that there is some evidence 
that the value of ae remains "constant for all laminates of all 
fiber reinforced/resin matrix composites...at least for what has 
been referred to as 'fiber of filament-dominated' laminates in 
glass/epoxy, boron/epoxy, and graphite/epoxy systems." There 
seems to be some difference, however, between tape (ae + 0.28" Еог 
AS/3501-5) and fabric (ao = 0.33" for HMF330C/34). The actual 
strength reduction is based on the ratio of the solid panel (PX@@) 
ultimate strength to that of each panel with the stress concentra- 
tion. The percent difference (> A) between calculated and actual 
strength reduction is [(calculated-actual)*(100/calculated)]. 
This value serves to compare the relative magnitude of observed 
strength among test specimens. A positive value of %А indicates 
a panel which demonstrated higher strength than predicted by the 
SCF computed by the LEFEA. Note the close correlation among FAPF, 
actual strength reduction and failure type. 
EL Types of Panel Failure 

There were two types of panel failure: (Type-1) delamina- 
tion at the point of highest stress concentration (0 - O° on the 
edge of the cutout) followed immediately by total failure and 
(Type-2) facesheet separation followed at some higher load by 


catastrophic failure.  Type-2 failures occurred far below the 
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expected stress level. A panel with a Type-2 failure not taken to 


to a complete failure was designated Type-2'. 


TABLE VIII 


COMPARISON OF EXPERIMENTAL RESULTS 


Panel | ----------------—- Ддҙ -------------------- Failure 
Desig Load (psi) SER Calc Actual % Type 
2×00 25,500 57,460 1.00 0.00 0.00 (RS 1 
PONY 17,000 30,000 99599 1.89 1.88 -0.5 Т! 
Jl l 18, 200 о 2 5 87 1.92 - 2.7 ЈЕ 
RR22 16,000 21,050 2. 39 187 2 7 -46.0 2 
БЕ! 1735809 2-29 2250 JL 51817] 2-05 -9.6 А 
RR42 L3 29500 21,900 2 1.55 2 -40.3 2 
RR 2, 3900 167000 27 1:37 ^2 а 2 
8511 17,000 317000 2.64 ей 85 abc Пі 
БӘСІ ТУ. = و‎ ٦ 20 (81 لس‎ +5.3 1 
RSSI 7,000 19,600 2.46 Isa, ** ШЕ; 2 
RH11 L9 00 29,960 Ваша je 1.92 —-1.6 1 
ЕН22 13,000 31,640 DI 0,۲ 1032 3.2 | 
RH E 2 218590 2752 esis 2267 -42.0 2 
RH42 2177000 36,990 2.64 1. ° 1.55 +17.1 l 
R451 lo, 500 317630 Zoe 1.88 Ia? +3 1 


Strength reduction calculation used as = 0.33 in. from LMSC test 


data. Panels marked with ** were not taken to total failure. 
Failure Types: 1 - Failure originates at Strain Concentration 
2 - Facesneet Separation & Buckle 
2' — Facesheet Separation & Buckle (not loaded 


to ultimate) 


It has been noted that in compression there exists a 
tensile Poisson stress (+0 _) which is greatest ater Ere 
highest stress concentration on the edge of of the cutout and 
tends to pull the plies apart. It was not possible from this 
experimental procedure to determine if ply delamination, inter- 
laminar shear stress or micro-mechanical fiber buckling was the 
initiator of the failure. In fact, failure probably resulted from 
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two or all three of these working together, possibly in conjunc- 
tion with both fiber and matrix flaws. Figure 5.5 shows a typical 
compressive panel (Tyoe-1) failure. All Type-1 failures were 
almost identical in appearance; they differed only in the ultimate 
load sustained. 

All Type-2 and Type-2' failures were also similar to eacn 
other; the facesheet began to pull away from the core at some 
point away from the edge of the cutout. This began with the 
formation of a small bulge or "bubble" which increased in total 
area and distance from the face of the core to inside surface of 
the facesheet. The initial separation was not visible until well 
into the load cycle, however, in some cases the FAPF may have well 
been the sound of the initial adhesive failure. Once the leading 
edge of the separation reached the cutout the panel failed 
Bpotally. 

Tyoe-2' failure was this facesheet separation not taken to 
total failure. The partially failed panels were removed from the 
test apparatus and subjected to non-destructive (NDI) and 
destructive inspections in an attempt to determine the possible 
cause of the core-facesheet separation. 

Panels RR22, RR42, RR51 RS51 and RH3l failed by facesheet 
separation. The stress-strain curves for these panels appear in 
the individual appendices and all clearly show the result of the 
facesheet separation--the slope of the curve dramatically in- 
creased.  Tnis was due to the decreasing panel stiffness and the 
picking-up of the load as the area of separation and facesheet 
curvature increased. 

When this failure mode appeared an additional FEA was 


considered necessary to examine in detail the core-facesheet 
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interface. A three-dimensional analysis was made of panel RR22 to 
determine if any significant out-of-plane stress was causing the 
separation. The results are given in Appendix O. The interface 
(the idealized adhesive surface) showed very low stresses in the 
+z, or out-of-plane, direction. From this it may be assumed that 
facesheet separation was the result of an incomplete or bonding 
process or other manufacturing error. 

The Type-2' failure panels were subjected to non-destruc- 
tive inspection using C-Scan and X-ray methods to attempt to 
locate the source of the defect(s). No obvious flaws or 
manufacturing errors were apparent. Panels virtually identical to 
those Type-2 and -2' failures were manufactured and tested under 
tne same experimental conditions. In each case the panels 
sustained Type-1 failure and carried an ultimate load into the 


-29 to -35 ksi range. 
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УТ. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

This study examined three geometric configurations of co-cured | 
reinforcement of graphite/epoxy honeycomb plates with circular 
cutouts subjected to uniaxial compressive loading and compared 
them to ldentically loaded unreinforced notched and solid plates. 
The test specimens were modeled using linear elastic finite 
element analysis (LEFEA) to analyze the strain field around the 
cutout. The objective of the study was to determine if a 
relatively simple, inexpensively manufactured reinforcement of a 
cutout could significantly reduce the stress concentration it 
induced, decrease the local strain and thereby increase the 
ultimate (failure) strength of the panel. 

Table IX is a summary of the important analytical and experi- 
mental results. The computed SCF is derived from the LEFEA. The 
predicted failure stress is based on the actual failure of the 
unreinforced panel (PO@%) and the analytical SCF. More than many 
any other experimental results, compressive failure in composite 
plates should be classed a stochastic function. It would take a 
number of identical panels of each configuration to arrive at a 
statistically significant predicted failure stress. However, 
from the data of this study, the average of the eight Type-1 
failures was 93.5$ of the predicted failure stress. The strip 
reinforcement (four Type-l failures) failed at 100.5% of the 
predicted applied stress. A three-dimensional linear finite 
element analysis of a typical Type-2 failure (Panel RR22) was 


attempted (see Appendix O). It failed, however, to predict the 
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actual failure loading or to provide a reason for the premature 


facesheet separation. 


The test program reported here confirmed tnat, properly used, 


the linear elastic finite element method provided an exceptionally 


accurate strain field representation even in a material with 


nonlinear response (see Appendices A-O). 


The failure stresses 


were harder to predict using linear methods, but this is hardly 


surprising considering the material is a composite and the loading 


ls compression. 


B. 


TABLE IX 


SUMMARY OF ANALYTICAL AND EXPERIMENTAL RESULTS 


Test Computed 


Panel 


ОСЕ 


Failure Stress (psi) 


Predicted 


Predicted Failure 
Load 


Type 
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Жаре 1“ а: сап ое seen that reinforcement 


reduced computed stress concentrations up to 20% in some configu- 


"۱۹. 


Тпе rexnftorcement added little more tnan l to 4$ 


additional weight to each panel. While it is difficult to 
directly compare the improvement reported nere with configurations 
іп an actual application in a large, complex structure, it is easy 
to see that for a small increase in weight a significant reduction 
іп stress concentration is both predicted and realized. 

Small amounts of graphite/epoxy reinforcing lamina(e) co-cured 
wlth thin composite sheets of the same material can significantly 
reduce stress concentrations and increase ultimate failure load. 
This reinforcement method involves some additional manufacturing 
effort, but it yields excellent strength-to-weight comparisons. 
The analytical results indicate that using several small 
reinforcement plies concentrated close to the cutout provides the 
most attractive strength-to-weight ratios. The strip 
configuration also gives excellent results and seemingly very 
oredictable failure levels. 

This experimental program reaffirmed the well-known fact that 
even minor manufacturing defects can be a severe problem in 
compression testing. Improper or incomplete bonding of the face- 
sheets to the honeycomb core can significantly affect the ultimate 
failure load in graphite/epoxy specimens. Іп five cases the 
facesheet began separating from the core at a point away fron ИЕ 
cutout. A "bubole" then formed reducing the facsheet's load 
resistance and transferring the load to the opposite, still intact 
facesheet. The panel then began to exhibit greatly decreased 
stiffness. As the load was increased, panel stiffness decreased 
in proportion--similar to Euler column buckling. 

It was not possible to locate the source of the bonding 
failure or even prove conclusively that improper bonding was the 


source of the premature failure. However, since the "bubble" 
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usually initiated at low load levels and at points well away from 
the stress concentration, bonding failure appears to be the most 
logical explanation. Prior to testing the panels were subjected 
to NDI which failed to discover any unbonded areas petween the 
core and facesheet. These failures could have been a case of weak 


or only partial adhesive bonding. 


C. RECOMMENDATIONS 

The research reported here investigated only a few of the 
possible reinforcement geometries. Any number of significant 
questions remain unanswered in the research reported here. 
Additional work is suggested in the following areas: 

Ши Еше-пег Testeng OF Reported Geometries 

Time and money limited testing to one specimen of each 
geometry. Several of the most promising reinforcement configura- 
tions (RR22, RR42, RH42, etc.) should be subjected to further 
testlng to obtain statistical confirmation of these results. 

The reaction of some of the strain gages remains 
unexplained, at least in part. For example, panel RS31 (Appendix 
H, Figure H.4) the gages closest to the cutout show points where 
an increase in load causes no corresponding strain increase. Ata 
higher applied stress the gage begins to react normally and 
stress-strain curve resumes an otfset but oarallel course (also 
see panel RH22, Figure K.4). 

2. Additional Reinforcement Geometries 

The three geometries reported here hardly exhaust the 
possibilities. Some additional promising configurations include 
oval (when the principal load direction is known or is predict- 
able), several different "wedding cake" methods and moving the 


Ew configuration closer to the cutout. 


POL 


3. The Effects of Reinforcement Stiffness 
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Reinforcement plies identical to the reinforced material 





was used in this study. It would be interesting to observe the 
effects of stiffer reinforcement such as laying G/Ep tape 
reinforcement O° to the applied load. 

4. Improvements to Experimental Methods 

A dense strain gage network next to the cutout on both 
sides of each facesheet may better explain the mechanics of 
failure.  Mucn closer load increments are necessary, lOO00 psi 
steps were not sufficient for a full explanation of the high 
strain notched panel response. 

A micro-photographic sequence of the stress concentration 
at the edge of the cutout at high load (starting atico To EOE 
might yield significant information on the way the graphite/epoxy 
panels fail in compression. 

The MTS macnine used in this research maintained a 
constant (or constantly increasing) load using an electronic feed- 
pacx loop. When the panel began yeilding, stiffness was reduced 
and the rate of nead travel increased in an attempet to maintain 
the indicated load. At failure, the head moved about 1/3" - 1/2" 
and crushed the panel. This precluded detailed examination of tne 
delamination at the stress concentration at failure. A constant 
displacement compression test machine is recommended in suosequence 


research. 
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APPENDIX À 


PANEL POØØ: ANALYTICAL AND EXPERIMENTAL DATA 


Panel 20۴٣) served as the basis for comparison between 
reinforced and unreinforced compression specimens. It had a cen- 
tered 1.00" diameter hole with no reinforcement around the cutout. 
Two identical specimens were produced, instrumented and tested. 
They both failed at the hole edge (Type-1) at an average apolied 
normal compressive stress (G) ОЁ -30,500 psi. 

The panel finite element model (mesn) is illustrated in Figure 
A.l. Figure 5.1 shows the distribution of strain around the 
cutout comparing open and closed-form computation methods. Table 
X gives tne (finite element) computed strain data around the 
cutout. Figure A.2 shows the correlation between computed strain 
(solid and dashed lines) and the experimentally measured strain 
(triangles) along the X axis at an applied normal stress of 
-10,000 psl. Figure 5.2 shows the correlation between open and 
closed-form analysis along the X axis. Table XI lists the 
computed values of the strain parallel to the applied load (Epbs-Y) 
and Poisson expansion (Eps-X) along the X axis. 

Note that in Figure A.2 between the 1.25 and 2.75 inch 
stations on the panel's X axis the indicated gage strain seems to 
alternate slightly up and down. Gages indicating nign were on the 
right side of the hole while those indicating slightly lower were 
on the left. In order to best illustrate the effect of the stress 
concentration, left and right side gages are superimposed on the 


right side of the cutout. It appears that the right side saw 
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about 1 to 2% higher strain than the left. The strain difference 
is attributed to either a very slight test fixture misalignment 
or a difference in panel length between each side of the hole 
amounting to about to somewhat less than 0.0005" 

Figure A.3 shows graphically the experimentally measured 
values of strain at different locations on the X axis. The 
numerical strain data are given in Table XII. The center of gage 
#l's resistive grid at x = 0.570" was 0.070" from the edge of the 
hole. The strain indicated was appropriate to the applied load 
taking into account the nonlinearities discussed in Chapter 5. 
Between -7,000 and -9,000 microstrain on gage +1, there appears to 
be a slight anomaly where the strain does not increase as fast as 
it nad up to that point, but it then appears to “catch 0٦ 
may be attributed to minor fiber failure, nonlinear load transfer 
or local delamination. This is a phenomenon that becomes much mGme 
apparent in the tests of reinforced panels. 

Figures A.4 through A.7 show the strain field contours at an 
apolied normal stress of -10,000 psi computed and plotted using 
DIAL. Figure A.4 is the full quarter panel which shows the effect 
of not loading the full width of the top edge. Figures A.5 
tnrough A.7 show the strains contours close to the cutout. The 
computed strains are very close to that for the ideal тип ин 
width panel. 

Figures A.8 through A.ll plot panel stress contours at the 
-10,000 psi loading. The type of element used in the LEFEA 
required the applied load to be input as a stress resultant (N, 
lb/in) and produced plots in the same units. Since all the 
reinforced panels had different thicknesses over their surfaces, 


the plots of stress resultant were not valid and are not given for 
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any other panels. Plots of stress are included to be compared 
with the classic notched plate solution to validate the analysis. 
Figure A.8 shows the full quarter panel with the stress 
concentration in the upper right corner due to the panel clamping 
modeling and the imposed boundary conditions. Figures A.8 and A.9 
show the maximum resulting stress parallel to the load. The 
maximum induced stress (at Û = 0°) is 31,100 psi which compares 
to 30,500 psi (Equation 2.3) for the finite-width panel. This is 
Just 3.6% over that predicted for an infinite plate. This minor 


difference is accounted for by the loading and boundary conditions. 
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Figure A.l Рапе1 РО@@: DIAL Finite Element Mesh. 
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TABLE X 


LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (—10,000 PSI). 
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Strain Comparison Along the X Axis. 
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Figure A.2 
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TABLE XI 


LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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Mircostrain vs. Cmpressive Stress. 
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Figure A.3 
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TABLE XII 
PANEL PO): SELECTED STRAIN GAGE VALUES DURING LOAD. 


Micrc-strain Indicated by gage: 
#1 


Load 
(PS1) 
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Panel POØØ: Eps-Y FEA Contours. 


Figure A.4 
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Panel POØ: Eps-Y FEA Contours Near the Cutout. 


Figure A.5 





Ерз-Х FEA Contours Near the Cutout. 
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Figure А.6 
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Panel POØØ: Eps-XY FEA Contours Near the Cutout. 


Figure A.7 
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APPENDIX B 
PANEL RR11: ANALYTICAL AND EXPERIMENTAL DATA 
Panel RRll was reinforced with one round co-cured ply of G/Ep 


concentric with the cutout placed on the outside of each facesheet. 


The reinforcement had the following dimensions: 


Shape: Round 
Inside Diameter: 1:00. in 
Outside Diameter: 2.24 ти 


Thickness (each): 0.014 in (1 ply) 
Area (each face): 3.16 in? 
Total Volume: 0.088 in? 


Net Cross Section: 0.035 in? 


The panel failed at the hole edge (Type-1) at an applied 
normal compressive stress of -29,950 psi. Based strictly on the 
failure of the unreinforced panel and the computed stress concen- 
tration factor of 2.65, the predicted failure was On = -34,500 
psi. Failure was at only 87% of the predicted load. There was no 
obvious reason for the early failure. No manufacturing errors 
were apparent on post-test visual or non-destructive inspection of 
the facesheet-honeycomb bonding. 

The finite element model (mesh) is shown in Figure В.1. The 
area of round reinforcement is denoted by the area inside the bold 
outline around the cutout. 

Figure B.2 compares the computed (finite element) strains 


around the cutout between the unreinforced panel (РО00) and RR11. 
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Table XIII gives the computed distribution of strains around the 
cutout in the Y and X directions as well as shear (Eps-X, Eps-Y 
and Eps-XY). 

Figure B.3 compares the finite element models (solid and 
dashed lines) and the experimental values (triangles) of strain at 
-10,000 psi applied normal stress. It shows the very close corre- 
lation between the analytically predicted and experimental strain 
with some ininor variation on either side of the panel. The LEFEA 
strain values are listed in Table XIV. The edge of the reinforce- 
ment extended to 1.12" in on the X axis. This is apparent from 
the figure as the inflection point in the direct compressive 
straia (solid line) where it abruptly begins to increase. 

Figure B.4 shows the stress-strain state during the load 
sequence from O to -30 ksi. Experimentally measured strain gage 
values are given in Table XV. At -16 ksi the gage next to the 
hole (x = -0.571") suddenly indicates a severe loss of local 
stiffness. This is reflected to a smaller, but no less dramatic 
degree in gage +3 on the other side of the cutout at x = 0.749". 
Gage #1 demonstrates what appears to be a continuous increase in 
local stirfness starting at -18 ksi; as the load increases the 
strain decreases. Between -22 and -24 ksi the strain is rapidly 
changing from compressive to tensile next to the hole. No visible 
buckling or delamination, which might help to explain part of this 
behavior, was noted next to the cutout under visual inspection. 

Figures B.5 through B.8 show the LEFEA computed strain 
contours at an applied normal compressive stress of -10,000 osi 
computed and plotted using DIAL. Figure B.5 is the full quarter 
panel with strain parallel to the applied load. This shows some 


strain contours at the top right of the panel illustrating the 


E 







effect of not loading the full width of the top edge. Pirgu а 
through B.8 (Eps-Y, Eps-X and Eps-XY) show the strains in detail 


close to the cutout. 
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Strain Comparison Around the Cutout. 
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TABLE XIII 


PANEL RRIl: LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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Strain Comparison Along the X Axis. 
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Figure B.3 
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TABLE XIV 


PANEL, RR11, LEFEA STRAIN DISTRIBUTION ALONG Til X AXIS (-10,000 PSI). 
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Microstrain vs. Compressive Stress. 


Panel RRll 


Figure B.4 
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Micrc-Strain Indicated by Gage: 
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SELECTED STRAIN GAGE VALUES DURING LOAD. 
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Panel RRll: Eps-Y FEA Contours. 
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Figure B.5 
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Panel RRl1: Eps-Y FEA Contours Near the Cutout. 


Figure B.6 
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Eps-X FEA Contours Near the Cutout. 


Panel RR11 


Figure B.7 
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Panel RRll: Eps-XY FEA Contours Near the Cutout. 


Figure B.8 


APPENDIX C 
PANEL RR22: ANALYTICAL AND EXPERIMENTAL DATA 
Panel RR22 was reinforced with two round co-cured plies of 


G/Ep concentric with the cutout on the outside of each facesheet. 


The reinforcement had the following dimensions: 


Shape: Round 
Inside Diameter: 1.00 in 
Outside Diameter: 2.24 iñ 


Thickness (each): 0.028 in (2 plies) 
Area (each face): 3.16 in? 
Total Volume: 0.176 in? 


Net Cross Section: 0.069 in? 


The panel failed by facesheet separation and buckling (Type-2) 
almost immediately upon initial application of the load. It 
failed totally at an applied normal compressive stress (0 ,) of 
-21,050 рві. Based strictly on the failure of the unreinforced 
panel and the computed stress concentration factor of 2.39, the 
predicted failure was On = -38,250 psi. This reinforcement 
configuration should have been among the most efficient: stacking 
the most additional thickness closest to the point of highest 
stress concentration. 

The finite element model (mesh) is shown in Figure C.l. The 
round area of reinforcement is outlined by the heavy lines next to 
the cutout. The reinforcement is two plies thick on the outside 


of both facesheets. 
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Flgure C.2 shows the comparison between finite element 
computed strains around the cutout between the unreinforced panel 
(РО00) and RR22 at -10,000 psi applied normal stress. Table XVI 
gives the computed values of the three strains (Eps-Y, Eps-X and 
Eps-XY) for the reinforced panel. 

Figure C.3 compares the analytical (solid and dashed lines) 
and experimentally measured strain values (triangles) at On = -10 
ksi. The alternating strain gage values between x = 1.0" and 2.5" 
indicate the small experimental difference between gages on 
opposite sides of the cutout. While the facesheet separation 
began at the onset of the load, it covered only a small area and 
the strain gages were on the opposite side of tne panel still 
giving reasonable indications at ОҒ = —-10 ksi. The edge of the 
reinforcement extended to x = 1.12". This can be seen clearly in 
Figure C.3 as the point where the strain along the X axis has an 
inflection point and begins increasing after a steady decrease 
moving away from the cutout edge. Table XVII gives the (finite 
element) computed distribution of strains around the cutout in the 
Y and X directions as well as shear (Eps-Y, Eps-X and Eps-XY). 

Figure C.4 shows tne stress-strain state during the load 
sequence from O to -20 ksi. Numerical strain gage data are given 
in Table XVIII. Gages #1, #3 and #4 all indicate a decreasing 
compressive strain rate with load application. It appears that a 
facesheet separated from the honeycomo core at or shortly after 
load application. When this is compared with the strain levels 
shown in Figure C.3, it appears that significant separation did 
not occur until after the -10 ksi load level. 

Begures C.5 through C-8 show the analytical strain contours at 


an applied normal compressive stress of -10,000 psi computed and 


p 


plotted using DIAL. Figure C.5 is the full quarter panel with 
strain parallel to the applied load. This shows some strain 
contours at the top right of the panel illustrating the effect of 
not loading the full width of the top edge. Figures C.5 through 
C.8 (Eps-Y, Eps-X and Eps-XY) show the strains in detail close to 
the cutout: 

Panel RR22 should have been the most efficient reinforcement 
configuration with the best ratio of volume-to-strength. While 
the volume of reinforcement was small, most of it was concentrated 


adjacent to the hole in the area of highest stress concentration. 


Note: After this research program showed premature panel 
failure due to facesheet separation, two additional RR22 panels 
were fabricated by LMSC using identical materials (a different 
lot, however) and methods. They failed at O, = -41.5 and -38.0 
xsi or an average 104% of the predicted applied normal stress of 


=33.3 кеш 
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REINFORCEMENT (2 plies/facesheet) 


DIAL Finite Element Mesh. 


Panel RR22 


Figure C.1 
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Strain Comparison Around the Cutout. 
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LEFEA STRAIN DISTRIBUTION AROUND 'PHE CUTOUT (-10,000 PSI). 
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Strain Comparison Along the X Axis. 
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PANEL RR22, LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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Microstrain vs. Compressive Stress. 


Panel RR22 


Figure C.4 
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TABLE XVIII 


SELECTED STRAIN GAGE VALUES DURING LOAD. 


PANEL RR22: 


#4 #15 


Micrc-Strain Indicated by Gage: 
# 3 


#1 
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Panel RR22: Eps-Y FEA Contours. 


Figure C.5 
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Panel RR22: Eps-Y FEA Contours Near the Cutout. 
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Figure C.6 
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Panel RR22: Eps-X FEA Contours Near the Cutout. 
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Panel RR22: Eps-XY FEA Contours Near the Cutout. 


NIIS 8 ٤ 
۱۱۱۱۱۱٢۹١ 


) ) ۹۵۳ 


ІП ЧИШ NOME 


] ٍا۱‎ 
IL. 1. ۷ 


Figure C.8 


=) 


147 


APPENDIX D 


PANEL RR31: ANALYTICAL AND EXPERIMENTAL DATA 


Panel RR3l was reinforced with one round co-cured ply of G/EP 
around the cutout on the outside of each facesheet. The rein- 


forcement had the following dimensions: 


Shape: Round 

Inside Diameter: 15:00 an 

Outside Diameter: 3-60 10 
Thickness (each): 0.014 in (1 ply) 
Area (each face): 9.39 in2 

Total Volume: 0.263 in? 


Net Cross Section: 0.073 in? 


The panel failed at the hole edge (Type-1) at an applied 


normal stress (Ø of -28,000 psi, only about 78% of that 


не, 
expected. Based strictly on the failure of the unreinforced panel 
and the computed stress concentration factor of 2.56 the predicted 
failure was 0, = -35,700 psi. 

The finite element model (mesh) is shown in Figure D.l. The 
area of round reinforcement is denoted by the area inside the bold 
outline around the cutout. 

Figure D.2 compares the analytical values of strain around the 
cutout from  0- 0° to 90? between the unreinforced panel (PO) 
and RR31 in the Y and X directions as well as shear (Eps-Y, Eps-X 
and Eps-XY). Table XIX lists these computed strains around the 


GUEGOUU: 
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Figure D.2 compares the finite element model (lines) and the 
experimentally measured strain data (triangles) of strain at 0, = 
-10,000 psi. It shows an almost perfect correlation between the 
analytically predicted and experimentally measured strain. The 
edge of the reinforcement extended to 1.80" in on the X axis. 

This is apparent from the figure as the slight inflection point 
where the strain begins increasing slightly. The computed strain 
field in an unreinforced panel (PO%%) is shown as dashed lines 
(Eps-Y & Eps-X). Table XX gives the values of the computed strain 
in the Y and X directions as well as shear. 

Figure D.4 shows the stress-strain relation during the load 
sequence from O to -28 ksi. Experimentally measured strain gage 
values are given in Table XXI. Other than a minor "glitch" at 
0, = ~21 ksi, no exceptional anomalies were noted. Gage #3 at 
x — -0.770" showed little increase in strain between -21 and -22 
ksi. This is not reflected in any of the other yage readings. 

Figures D.5 through D.8 show the strain contours at an apolied 
normal stress of -10,000 psi computed and plotted using DIAL. 
Figure D.5 is the full quarter panel with strain parallel to tne 
۲۰٠٦۰۰۹ 160. As before, the strain contours at the top right of 
tne panel are due to the effect of not loading the full width of 
the top edge. Figures D.o through D.3 (Eps-Y, Eps-X and Eps-XY) 


show the strains in detail close to the cutout. 
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REINFORCEMENT (1 ply/facesheet) 


Figure D.l Panel RR31: DIAL Finite Element Mesh. 
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Strain Comparison Around the Cutout. 


Panel RR31 


Figure D.2 
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‘TABLE XIX 


LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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Strain Comparison Along the X Axis. 


Panel RR3l1l 


Figure D.3 


TABLE XX 


LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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Microstrain vs. Compressive Stress. 


Panel RR31 
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TABLE XXI 


FED STRAIN GAGE VALUES DURING LOAD. 
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Panel RR31: Eps-Y FEA Contours. 


Figure D.5 
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Panel RR31: Eps-Y FEA Contours Near the Cutout. 


Figure D.6 





158 


5 


0+ + 7 
А ie 


Ee 
СОЗО 


а 100 
а 100 
ии 
ШЕШІ 
и 1001 

МІН 


| [outs 
JON Vel 


[ ши, 





г. 
ОИ t ООО а | 








— — — — ee 


کس ———————————————— 


"n - XUW 
O NIN 


П а 
| `Ú / 


го IDORO й 


4701 5 
ho || 
ИВ ( 
(ра C 
СИ ШЕ | 
cx ۲۷ 

ННОТМО 1 


IN IW 133 
ЧООТМИ ) 


IN WIT I 
ша ШЫ! 





Panel RR31: Eps-X FEA Contours Near the Cutout. 


Figure D.7 
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Figure D.8 





APPENDIX E 
PANEL RR42: ANALYTICAL AND EXPERIMENTAL DATA 
Panel RR42 was reinforced with one co-cured ply of G/EP 


Esuncentric with the cutout on tne outside of each facesheet. The 


reinforcement had the following dimensions: 


Shape: Round 
Inside Diameter: OO AR 
Outside Diameter: ЗОО un 


Thickness (each): O OPO ARE DIY) 
Area (each face): 6.28 in? 
Total Volume: 0.176 in? 


Net Cross Section: 0.112 ia^ 


The panel failed by facesheet seoaration and buckling (Type-2) 
at an applied normal stress (0,) of -13,500 psi, less that 34% of 
the expected value. Based strictly on the failure of the unrein- 
forced oanel and the computed stress concentration factor of 2.29 
for RR42, the predicted panel failure маз 0, - -40,000 psi. 

The finite element model (mesh) is shown in Figure E.l. The 
area of round reinforcement is denoted by the area inside the роја 
outline around the cutout. 

Figure E.2 compares the LEFEA values of strain around the 
cutout between the unreinforced panel (POUJ) and RR42 in the Y and 
X directions as well as shear (Eps-Y, Eps-X and Eos-XY). Table 


XXII lists these computed strains around the cutout. 


Flgure E.3 compares the finite element model (lines) and the 
experimentally measured strain data (triangles) of strain at ба = 
—10,000 psi. It shows very poor correlation petween the 
analytically oredicted and measured strain. At least one face- 
sheet separated from the honeycomb core over a significant area 
prior to the -l0 ksi test point. The strain field at that applied 
stress was little more than 80% of that predicted. The computed | 
strain field in an unreinforced panel (PO@@) is shown as dotted | 
lines (Eps-Y & Eps-X). Table XXIII gives the values of the LEFEA 
computed strain in the Y and X directions as well as in shear. 

Figure E.4 snows the stress-strain relation during the load 
sequence from 0, = 0 to -21.9 ksi (failure). Measured за 
gage values are given in Table XXIV. The strain, particularly in 
gages #3, #8 and #14, indicate that the panel stiffness decreased 
from the initiation of the load. The notable difference in slope 
between gages #3, #8 and #14 and gage #1 next to the cutout 
indicate that the separation "bubble" occurred амау Frommer. 
cutout. Failure seemed to occur when the bubble's edge reached 
26761106611٤ 

This type of failure probably indicates that at least one 
tacesheet was improperly bonded to the honeycomb core.  Non- 
destructive inspection before and after testing did not indicate 
unbonded areas between the facesheet and the core. I suspect, in 
this and other panels that failed in a similar manner, that tne 
adhesive was in place but either weak from an improper mixing or 
aging or was applied too thinly. 

Figures £.5 through E.8 show the strain contours 88ب ج‎ ыа 
normal stress of -10,000 psi computed and plotted using DIAL. 


Figure E.5 is the full quarter panel with strain parallel to tne 
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applied load. As before, the strain contours at the top right of 
the panel are due to the effect of not loading the full width of 
the top edge. Figures 6.6 through E.8 (EPS-Y, EPS-X and EPS-XY) 


show the strains in detail close to the cutout. 
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TABLE XXII 


ГЕРБА STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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Figure E.3 
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TABLE XXIII 


PANEL RR42: LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (—10,000 PSI). 
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Microstrain vs. Compressive Stress. 


Panel RR42 


Figure E.4 
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Panel RR42: Eps-Y FEA Contours. 


Figure E.5 
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Panel RR42: Eps-Y FEA Contours Near the Cutout. 


WHS ХНИ | 

jJ o£ '0--NIMW 

' 1000 °0 / 

БЕЛІП» 9 
ОТО P 
ше И H 
1007 °0 f 
IAG? OO С 
(ООЕ ТО | | 


NUTS ۰۹ 
ИПОТНО ) 


) 71 
рат ипи1мо ) 


Figure E.6 


E T6 IN1N111 
۶+" " ۶“ ۹ 


172 


— 


-ad 


~ 
— 


ant 
1 
4 
4 


гу 
3H 


— . — > —— — — - — — 3ц_——_-ы— 


, 


C1 
~ 


- 


Мн ы ا‎ ees — —— - 
~ 
А 
. 
— 





мВ li = 
(кә ut 


rr. 
0)0 ئ٦‎ 


0710 


_ 


. 
— 


РИ Хе | 





КЕЕ ЧП ШЕ 
ir Svp ПИЯ 


ШІ 

ОТО (1 
11106 f J 
ane tv) Я 


| 
| 
ИИИ | | 
IH °! í) 
0101.0 h 
ОНО 72 
| 
| 


100: 700 / 
oDe TD 4 
ПТО ٦ 
» 1700070 " 
ШЕГІ? f 
NG p 
iine coo l 
MITES X HINT 
HOOmLNOO 


IGR N IH lI 
E I 


Е, ПМ МЕТ | 
ШЕГІ 


Panel RR42: Eps-X FEA Contours Near the Cutout. 
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Figure E.7 
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Panel RR42 


Figure E.8 


APPENDIX T 
PANEL RR51: ANALYTICAL AND EXPERIMENTAL DATA 
Panel RR51 was reinforced with one round co-cured ply of G/Ep 


in a concentric with the cutout placed on the outside of each 


facesheets. The reinforcement had the following dimensions: 


Snape: Round 

Inside Diameter: 111 ٤ 

Outslde Diameter: 4.60 іп 
Thickness (each): О ан дот Бед 
Area (each face): 15.83 in? 

Total Volume: 0.443 in? 

Nem cross section: 0.101 in? 


The panel failed by facesheet separation and buckling (Type- 
2') and was taken only to an applied normal stress of -15,000 psi, 
not to total failure. The series of premature panel failures due 
to facesheet separation required an intact panel for testing. 
Subsequent non-destructive testing showed the separation, but 
could not determine the reason for it. It is suspected that the 
adhesive, while properly applied, was not properly mixed or was 
overage. Based strictly on the failure of tne unreiaforced panel 
ema the computed stress concentration factor of 2.51, failure was 
predicted at about 0, = -36,400 psi. 

The finite element model (mesh) is shown in Figure F.l. The 
15.83 square inch area of the one-ply reinforcement is outlined by 


E Dold lines next to the cutout. 


i5 


Flgure F.2 compares the three (finite element) computed 
strains around the cutout between the unreinforced panel (PO) 
and RR51. These computed strain values are listed in Table ۰ 
Note the very significant decrease in the strain due to the rein- 
forcement at the point of highest stress concentration ( 60 = 0?) 
compared to the unreinforced panel. A significant decrease in all 
three strains can be seen around the hole from 0 to 90 degrees. 

Figure F.3 compares the LEFEA computed (solid and dashed lines) 
and experimental strains (triangles) in the Y and X (poisson 
expansion) directions in the panel and shows that there was a 
great disparity between opposite sides of the nole on the same 
facesheet. One side showed much higher strain than predicted at 
ga = -10,000 psi. This is due to load transfer trom TE 
With the buckled facesheet. The edge of the reinforcement can be 
seen in the figure by the very slignt inflection point at x = 
2.3". The effects of the one-ply reinforcement iS apparent in the 
far-field as a significant decrease in computed Eps-Y compared to 
the unreiaforced panel (POØØ). Table XXVI gives the computed 
values of the strains along the X axis. 

Figure F.4 snows the stress-strain state during the load 
sequence from 0 to -16 ksi. Experimentally measured strain gage 
values are given in Table XXVII. Gage #1 апа #4, on either side 
of the hole, show a positive slope of the derivative AQ/AE . 
This is unusual and indicates some panel-honeycomb seoaration 
close to the hole. 

Figures F.5 through F.8 show the strain contours at an om = 
-10,000 psi computed and plotted using DIAL. Figure F.5 is the 


full quarter panel with strain (Eos-Y) parallel to the applied 


ВИ 





load. Figures F.G through F.8 (Eps-Y, Eps-X and Eps-XY) show the 
strains in detall close to the cutout. 

The result of the wide reinforcement is to effect a thicker 
overall panel. A separate LEFEA of a [02,+45,90,соге]. panel 
without reinforcement had an almost identical strain field near 


the cutout. 
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Figure ١1 Panel RR5l: DIAL Finite Element Mesh. 
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Strain Comparison Around the Cutout. 


Panel RR51 


Figure F.2 
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PANEL RR51: LEEFEA SPRAIN DISTRIBUTION AROUND THE CUTOUT (—-10,000 PSI). 
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Strain Comparison Along the X Axis. 


Panel RR51 


Figure F.3 
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PANEL RR51: LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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Microstrain vs. Compressive Stress. 


Panel RR51 


Figure F.4 
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Panel RR51: Eps-Y FEA Contours. 
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Figure F.5 
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Figure F.7 
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Eps-XY FEA Contours Near the Cutout. 
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Panel RR51 


Figure F.8 


APPENDIX G 
PANEL RS11: ANALYTICAL AND EXPERIMENTAL DATA 
Panel 2511 was reinforced with one square co-cured ply of G/Ep 


Ewcentric with the cutout on the outside of each facesheet. The 


reinforcement had the following dimensions: 


Shape: Square 
Inside Diameter: IOS E LAO 
Length & Width: 2990 3m 


Thickness (each): ОА (I ply) 


Area (each face): 23020 in? 
Total Volume: 0.088 in? 
2 


Net Cross Section: 0.028 ја 


The panel failed at the hole edge (Type-1) at an applied 
Ш Па! stress of -31,000 psi, about 90% of the load predicted. 
Based strictly on the failure of the unreinforced panel and the 
computed stress concentration factor of 2.64, tne predicted 


failure was O0, = -34,600 psi. 


n 
The finite element model (mesh) is shown ian Figure G.l. The 
Square area of reinforcement is denoted by the heavy outline 
па the cutout. 
Figure G.2 compares the three (finite element) computed 
strains around the cutout between the unreinforced panel (POØØ) 


and RS11. These computed strain values are listed in Table 


XXVIII. 


09 


Figure G.3 compares the computed (solid and dashed lines) and 
experimentally measured (triangles) strain along the X axis in the 
Y and X (poisson expansion) directions in the panel and shows the 
excellent correlation between analytical and experimental strain 
at 0, = -10,000 psi. There was some minor strain variation 
between the left and right sides of the hole. Both are 
represented in the figure as strain gage values on the right side. 
The difference was small, but visible. The outside edge of the 
reinforcement can be seen in the figure as an inflection point in 
the direct compressive strain where it begins increasing at x = 
1.1". Table XXVIX gives the computed values of the strain along 
the X axis. 

Figure G.4 is the stress-strain state during the load sequence 
from gt = 0 to -30 ksi. Experimentally measured strain values 
are given in Table XXX. Up to about ~20 ksi the gage next to the 
hole (x = +0.561") shows an almost linear stress-strain relation. 
Gage #4 at x = -0.737" snows the expected degree of loss in local 
stiffness up to -20 ksi (see Table VI, Figure 5.4 and section III 
A.l for a discussion). Gage #1 indicates, starting at about -20 
ksi, what at first appears to be a slow but continuous increase in 
local stiffness indicated by a decreasing strain rate. At 
corresponding stress values gage $4 indicates an increasing s مه‎ 
rate. These gages were on opposite sides of the panel. What is 
actually happening is tnat 0,--28 көбі the right side ٣ 


panel is showing Significant matrix degradation and the load is 


being transferred to the left side of the panel next to the cutout 


Gage xl is probably not indicating the true state of strain Under 


its grid. Since the panel was evenly loaded and well constrainei 


0 


- а и аа — 


the strain at each side of the cutout 90° to the applied load 
should have been almost identical up to failure. 

Figures G.5 through G.8 show the strain contours at Qu = 
-10,000 psi computed and plotted using DIAL. Figure G.5 is the 
full quarter panel with strain (Eps-Y) parallel to the applied 
load. Figures G.6 through G.8 (Eps-Y, Eps-X and Eps-XY) show tne 


Strains in detail close to the cutout. 


1 
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Figure G.1 Panel RS11: DIAL Finite Element Mesh. 
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Figure G.2 
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TABLE XXIX 


PANEL RSl11, LEFEA STRAIN DISTRIBUTION ALONG ‘THE X AXIS (—10,000 PSI). 
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Figure G.4 
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TABLE XXX 


SELECTED STRAIN GAGE VALUES DURING LOAD. 


PANEL RS11 
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Panel RS11: Eps-Y FEA Contours. 
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Figure G.6 
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Panel RSll: Eps-X FEA Contours Near the Cutout. 
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Panel 2511 


Figure G.8 


APPENDIX H 
PANEL RS31: ANALYTICAL AND EXPERIMENTAL DATA 
Panel RS31 was reinforced with one square co-cured ply of G/Ep 


concentric with the cutout on the outside of each facesheet. The 


reinforcement had the following dimensions: 


Shape: Square 

Inside Diameter: +77 

Length & Width: 30 тп, 
Thickness (each): Dog при ply) 
Area (each face): SES in? 

Total Volume: 0.265 in? 


Мес Gross Section: 707062 in 


The panel failed at the hole edge (Type-1) at an applied 
normal stress of -32,550 psi, 89$ of the ultimate load predicted. 
Based strictly on the failure of the unreinforced panel and the 
milled stress concentration factor of 2.51, failure was 
oredicted at 0, 7(-36,400 psi. 

The finite element model (mesh) is shown in Figure Н.1. Тһе 
square area of reinforcement is denoted bv the bold lines next to 
pae cutout. 

Figure H.2 compares the three (finite element) computed 
Strains around the cutout between tne unreinforced panel (PO@¢G) 
and RS31. These computed strain values are listed in Table XXXI. 

Figure H.3 compares the computed and experimental strains in 


tne Y and X (poisson expansion) directions in the panel and shows 
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an almost perfect correlation between analytical and experimental 
strain at go = -10,000 psi. There was virtually no strain varia- 
tion between the left and right sides of the hole. The edge of 
the reinforcement is somewhat difficult to see in the Figure H.3 
as a slight inflection point at about x = 1.6". Table XXXII gives 
the LEFEA computed strain values along the X axis. 

Figure H.4 graphically shows the stress-strain state during 
the load sequence from 0, = 0 to -32 ksi. The experimentally 
measured strain gage values are given in Table XXXIII. Up to 
about 20 ksi gage #1 next to the hole (x = +0.553") shows an 
almost linear stress-strain relation. At or just above -20 xsi, 
however, there appears wnat seems to be a sudden decrease in 
strain rate on the right side of the cutout which just as suddenly 
ends at -23 ksi where the previous stress-strain ratio resumes. 

I believe that this is, instead, a transfer of very localized 
stress (or the load oath) away from the area next to the cutout to 
some other path in the field or possibly the opposite facesheet. 

t 1s important to note that gage #2 on the left side of the 
cutout at x = -0.597" shows no corresponding increase in strain 
that would oe caused by the transfer of load. Gage #3 at x = 
-0.660" shows some correspondence with gage 31 degree of loss in 
local stiffness up to -20 ksi. If it were a malfunction of the 
strain gage or a oartial debonding from the surface of tne 
composite the strain rate would change. 

Another anomaly occurs at p above -27 xsi where the "stair- 
step" phenomenon occurs again. Gage #3 at х = 0.606" reflects 
what is occurring next to the cutout edge. This 15 not ru K 
gages #6 and zll at x 2 -1.48" and 2.54" respectively; they 


reflect only the expected stress-strain relation. At -30 ksi gage 
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il indicated a rapidly increasing rate of strain and subsequently 
failed. 

Figures H.5 through H.8 show tne strain contours at On = 
-10,000 psi computed and plotted using DIAL. Figure H.5 is the 
full quarter panel with strain (Eps-Y) parallel to the apolied 
load. Figures H.6 through H.8 (Eps-Y, Eps-X and Eps-XY) show the 


Eamus in detail close to the cutout. 


PANEL R331: SQURRE REINFORCEMENT 
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Figure H.l Panel RS31: DIAL Finite Element Mesh. 
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TABLE XXXI 


LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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TABLE XXXII 


PANEL RS31: LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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TABLE XXXIII 


PANEL RS31: SELECTED STRAIN GAGE VALUES DURING LOAD. 
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Figure H.5 
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Panel RS31: Eps-Y FEA Contours Near the Cutout. 
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Eps-X FEA Contours Near the Cutout. 
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Eps-XY FEA Contours Near the Cutout. 
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APPENDIX 1 
PANEL RS51: ANALYTICAL AND EXPERIMENTAL DATA 
Panel RS51 was reinforced with one square co-cured ply of G/Ep 


concentric with the cutout on the outside of each facesheet. The 


reinforcement had the following dimensions: 


Shape: square 
Inside Diameter: lOO. in 
Length & Width: ДИК тя 


Thickness (each): 0.۰014 In (1I ply) 
Area (each face): 10.025 in 


Total Volume: 0.449 in? 


No 


Net Cross Section: 02087 in 


The panel failed at the hole edge (Type-1) at an applied 
normal stress of about -16,000 psi. Based strictly on the failure 
of the unreinforced panel and the computed stress concentration 
factor ot 2.46, however, the failure should have been close а = 
237,000 psi. 

The finite element model (mesn) is shown in Figure I.l. The 
square area of reinforcement is denoted by the heavy outline 
around the cutout. 

Figure I.2 compares tne three (finite element) computed 
strains around the cutout between the unreinforced panel (PO) 
and RS51. These computed strain values are listed in Table XXXIV. 

Figure I.3 compares the computed (solid and dashed lines) and 


experimental (triangles) strains in the Y and X (poisson 


Du 


expansion) directions in the panel and shows excellent correlation 
petween analytical and experimental strain at -10,000 psi applied 
normal stress. There was some minor strain variation between the 
left and right sides of the hole. The edge of the reinforcement 
is very difficult to see in the figure as a very slight inflection 
point at about x - 2.0". Table XXXV gives the analytical strain 
values along the X axis. 

Figure I.4 graphically snows the stress-strain state during 
the load sequence from O, = 0 to -21 ksi. Experimentally 
measured strain gage values are given in Table XXXVI. Up to about 
-16 ksi gage #l next to the hole (x = 40.553") shows the ехрес+ей 
almost linear stress-strain relation. However, just above ga = 
-16 ksi up to -19 ksi there begins a apparent loss in stiffness on 
the right side of tne cutout which suddenly ends at -19 ksi wnere 
the strain next to the cutout drops to almost zero. Note that 
даде #3 on the right side of the cutout at x = 0.666" snows no 
corresponding increase in strain that would be caused by an 
increase in local stress near the cutout due to a shift in tne 
load path. Gages #5 and #ll at x = 1.44" and -2.76" respectively 
reflect only the expected stress-strain relation. Тпіз сап be 
explained by gages #1 and 3 showing the effect of a gage under 
compression when the facesheet under it suddenly buckles outward. 
The result was a near zero strain indication. It is q1£Ficuim 
see, but there is an appreciable increase in strain rate indicated 
in gages #5 and #11 at O, = -20 ksi. This confirms that ЕВЕ 
a sudden increase in load in an area relatively far from the 
cutout, just as would be expected when the material close to the 
cutout begins to fail and the load paths are displaced away from 


it increasing the stain in the far-field. 


Figures I.5 through I.8 show the strain contours at O 


-10,000 psi computed and plotted using DIAL. Figure I.5 is the 


full quarter panel with strain (Eps-Y) parallel to the applied 


load. Figures I.6 through I.8 (Eps-Y, Eps-X and Eps-XY) show the 


strains in detail close to the cutout. 
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Figure I.l Panel RS51: DIAL Finite Element Mesh. 
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Strain Comparison Around the Cutout. 
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TABLE XXXIV 


PANEL RS51: LEPEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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TABLE XXXV 


LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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SELECTED STRAIN GAGE VALUES DURING LOAD. 
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Figure I.5 
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Figure I.8 


APPENDIX J 


PANEL RH11: ANALYTICAL AND EXPERIMENTAL DATA 


Panel RH11 was reinforced with one co-cured ply of G/Ep in the 
shape of two strips on either side of the cutout on the outside of 
each facesheet offset 0.50 inch from the edge of the cutout. The 


reinforcement had the following dimensions: 


Shape: S rio 
Length: 0 ۷ 
Width: 7 


Thickness (each): ШІ па (1 olv) 
Area (each Ғасе): 3.14 in? 
Total Volume: 0.088 in? 


Net Cross Section: 0:056 in? 


The panel failed at the hole edye (Type-1) at an applied 
normal stress (9,) ОЕ -29,960 psi. Based strictly on the failure 
of the unreinforced panel and the computed stress concentration 
factor of 3.02 (which was very slightly higher than the unrein- 
forced panel), failure was predicted at б, = ی۶‎ The 
panel failed witnin O.6$ of the oredicted ultimate load. 

The finite element model (mesh) is shown in Figure J.l. Tae 
area of the strip reinforcement is outlined by the heavy 11125 
offset 0.50 inch to the right of the cutout edge. 

Figure 2.2 compares the three (finite element) computed 
strains around the cutout between the unreinforced оапе1 (2089) 


and RHll. These computed strain values are listed in Table 


Zol 


XXXVII. There is no significant decrease in the strain qü отаца 
reinforcement. A very slight increase may be seen in Eps-Y near 
the O degree position (on the X axis). This increase in strain is 
due to the shifting of load paths to either side of the reinforce- 
ment. Тһе slight load path shift toward the cutout acted to 
Slightly increase the SCF. 

Figure J.3 compares the computed and experimental strains in 
the X and Y (poisson) directions in the panel and shows almost | 
perfect correlation between analytical and experimental strain at 
С = -10,000 psi. There was virtually no measured strain 
variation between the left and right side of tne hole. The edge 
of the reinforcement can not be seen in the figure; the effects of 
reinforcement 1S a decrease around the reinforcement and a subtle 
increase near the hole and in the far field (where x = 2.50") 
compared to the unreinforced panel (РОДФ). Taole XXXVIII gives 
the computed values of the strains along the X axis. 

Figure J.4 shows the stress-strain state during the load 
sequence from 0 to -29 ksi. Experimentally measured strain values 
are given in Table ХХХІХ. Up to -23 ksi all gages indicated a 
normal stress-strain state. From -23 ksi there was a dramatic 
change; first gages #1 and +2 showed a load transfer from the area 
next to the right edge of the cutout to the left side. Then 
suddenly tne roles were reversed and gage 72 (х = -0.583") 
indicated a load transfer to the right side of the cutout. Gage 
il (x = +0.568") snows a tremendous strain increase, o£ ane зса1= 
on Figure J.4, as high as 12,800 microstrain (Tables a The 
effect of the load transfer is apparent on gage #6 (x = +0.681); 
it reflects the increased Load away from the cutout edge on the 


right side of the hole. It apoears that the fibers on the right 
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side of the cutout began to buckle on the micro-mechanical level 
very near the edge. 

Figures J.5 through J.8 show the strain contours at 6 = 
-10,000 psi computed and plotted using DIAL. Figure J.5 is the 
full quarter oanel with strain (Eps-Y) parallel to the applied 
load. Figures J.6 through J.8 (Eps-Y, Eps-X and Eps-XY) show the 


strains in detail close to the cutout. 
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Strain Comparison Around the Cutout. 


Panel RHll 


Figure J.2 


TABLE XXXVII 


LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 


PANEL RHlL 


ЕР = EP Sak 


۴۴ 


СООКО 


Y 


COORD 


X 


МОРЕ 


77 СП СО СО СУ СОСО ДГ СОСЈОСОЈСОЈСКИМОМ СТ ٭ے‎ 
OOOO OCD O 5 = > = = 
| 1 | 1 I | J в IS 
53 22303 272 2] 213 21 23 #2 1:3 ہی ۰۹ تج‎ (1 AME E O O 
N= MO тусо оу ус Ре с Се Со Оре Sea NO 
Q Ye C2OCNLO 00 CO CN e CÓ ZPO ([— OY یی‎ CA CY O CN 
~M OnE NNS p Ns тото 
= Cry LO IOI یہی اج سم سم چہ‎ rry cr Cy cN CN — CD es 
e 8686 • $9 $ % $ $ о % 8 8 % ç ٭ 9 ?9 @ ق‎ 9 9 
СОЗ ОСУ ВУ ВЕЗУ Озе = = 

кк, ү, и 


cN NCNCN CN 070) Zr 0 0n 0n) ^n سم وب چ- سب‎ r y CN CN СУ 
Ооо сз О) сува > > = = تہ ری دہ زی‎ 
Il (I l| | ٘ ۱ً 1 7707 1+ +-+ 1:1:۵۲ 
£1 Fx Bd O O CG O FO FO BL ESL Fa p E ES C3 £2 E£ EC] 27121 
Гросо го те N= О ПИ О 9 22 С0 О о > 
CC ON TEINS GCG TOONS O TIO SONI о 
جح‎ C gD CED ہے‎ CO اوج جن تج‎ eS Gees | ele 
стз тт» =” «== [=ч (СУС Уе Г) ТТС) з се с: [хОу чу пе ти 
t€ ç ° ° 9 $ % ос о о о с с о о об с ос о ¢ ¢ 
ОС СЗО Савез среза ӘБ 
к ү |с! 


MOYO) OVS Se SOMA CIN GG پچ‎ 
Се оО ج‎ 
гурт Г 7٦ 
ВА ЕУ 2 ء١ د6‎ i) е (صع)‎ Гей кеф Со е рез Сет Ге Се بت‎ ۹ 
LO =F ins I= NFS BONO Sia = oS = 
заПа 5 uw ADDY FE ема ал өТө ее, лел Кә, 
دہ بے یج یہ عد یت ف‎ SIN رت ح نہ کر یہی‎ СЧ 
COLI SS COON ON) SONS OO Е ЧОЧО S 
eee °° 8 @ h 9 $9 о о с 9 % % е 9" t9 $9 © ¢ 
تے‎ c C e = e ٹک ری ری‎ 
РЕГ Рини Пи 


=> суд > S S ہے اج کے بے‎ те 
— C СО  ت کی بے بے بے دی تی‎ сос > css 
+ + + = + + + + + + + + ^ لب‎ + + | | + 
AOA AIA جج‎ 621 (14 212417] 6174... 36374 541 
СУ СО Со ОПИСИ СКОР СЛС СВ abe CISCO 
оосо) ео Р Осо: Mene TOCNO 
ОСОО ко CN O OL осел ری‎ O C O 
(ус? гу Оо — s => 
9 9 • • • • г % • 9 9 11 в v ° е • 


O ¢$ Ff è 
CCC CC O S = = = = 


зоо осоо ہج ے چ بے 2ے ہ ے ہے‎ 22 2 
Осе со о а COE CE 
+ ++ + ++ + ++ + + + + + + + + + + 
Ed EE بحا لع لعا لعا لعا‎ £n ыы ыы حر رع‎ E C3 E83 E Ee] £a 
DAA Го ле <> cN ONG осал У ЕСО 
ОСЕ" c = суш 2 22 


DAD п СУЧ СиО о СЧ ио ~ ОССО 
сг sJ с 


ооо 


C c So O S ООО ОСО СС е е =] 


— CN) cr) — tr) O I=* O GX OD ہے‎ чес Сосо کی‎ 


236 


(S9u9ut) sixy-x IƏued 








(peo10jure1ug) uieJ1 S 1vig 
(pe210jureg) ure ns туја 


шмәшәіпв геру әбес) 816 y 


۲8٢۷۲٦888 +۶٤5 





Адешш/с 40 sixy ہ6]10071010۱]‎ 6001۰7 ۱٥٢١۱6-1۷ 
(^S-) 559. 45 ом5зојашој 1564 000۶0٥۵6٥٢٥٥ د8ل‎ 
INOIND Jena раозојшен ас (инн |9эчез 


ISd 0000- 





000} 


0006 


ООО? 


(ui/ui) uie S -0191 


Strain Comparison Along the X Axis. 
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Figure J.3 
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TABLE XXXVIII 


PANEL RHl1l, LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (210,000 PSI). 
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Figure J.4 
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Figure J.6 
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Eps-X FEA Contours Near the Cutout. 
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Figure J.7 
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Figure J.8 


APPENDIX K 


PANEL RH22: ANALYTICAL AND EXPERIMENTAL DATA 


Panel RH22 was reinforced with two co-cured plies of G/Ep in 
the shape of two strips on either side of the cutout on the 
outside of each facesheet offset 0.50 inch from the edge of the 


cutout. The reinforcement had the following dimensions: 


Shape: SIT» 
Length: io? “ale 
Width: ۹۲ 


Thickness (each): 0255 (2 ply) 


Area (eacn face): ЭА іп2 


WW 


Total Volume: 0۰٦76 in 


Шер Gros: Section: 0.112 in? 


Tne panel failed at the hole edge (Type-1) at an apolied 
normal stress of -31,460 psi. Based strictly on the failure of 
the unreinforced panel and the computed stress concentration 
factor of 2.91, the predicted failure was б. = –31,500 р51. The 
actual failure was within 0.4% of the pcedicted ultimate load. 

The finite element model (mesh) is shown in Figure К.1. The 
area of the strip reinforcement is outlined by the heavy lines 
ге“ 0.50" to the right of the cutout's edge. 

Figure K.2 compares the three (finite element) computed 
strains around the cutout between the unreinforced panel (PO) 


and RH22. Tnese computed strain values are listed in Table XXXx. 


There is only a small decrease in the strain due to the 
reinforcement. 

Figure K.3 compares the computed (solid and dashed lines) and 
experimentally measured (triangles) strains in the X and Y 
(poisson) directions in the panel and shows an excellent 
correlation between analytical and experimental strain at -10,000 
psi applied stress. There was some very slight measured strain 
variation between the left and right side of the hole. The exact 
edge of the reinforcement can not be seen in the figure. The 
effects is a significant strain decrease (compared to the unrein- 
forced panel, PO) under the reinforcement, a small decrease near 
the hole and a slight increase in the far field (where x - 

2.50"). Table XLI gives the computed values of the strains along 
the X axis. 

Figure K.4 snows tne stress-strain state during the load 
sequence from ша = 0 to -30 Ksi. Experimentally measured strain 
gage data are given in Table XLII. Up to about -9 ksi all gayes 
indicated a normal stress-strain state. At -9 ksi gage #l (x = 
-0.570") began showing decreasing reaction to tne applied load. 
At -18 Ksi, whicn coincided with the first audible ply failure 
(FAPF), the strain at gaye #1 showed virtually no cnange uo to -22 
ksi. At -23 ksi, however, the strain suddenly doubles from 3100 
=o 6500u € and resumes its normal stress-strain ratio. At -9 ksi 
it appears that the load path is being diverted away from the 
right side of the cutout to the left. Gage #2 (x = -0.614") 
Aemonstrates an increased strain rate from -9 to -25 ksi when it 
increases significantly. Gage +1 failed above -27 ksi at about 
10,0004 € while Gage #2 continued to give reliable output uo to 


almost 12,000u&. It can ve assumed from the response of ygages ia 
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٤24,1‏ ٹیپ 


the far-fleld that the stress-strain response of tne panel аз а 
whole remained constant with a slight decrease in stiffness with 
increased loads. The response of the facesheets close to tne 
cutout show a very different response. It appears that there is a 
Significant transfer of load from one side of the cutout to the 
otner and to the opposing facesheet. 

Figures K.5 through K.8 show tne strain contours at а = 
-10,000 psi computed and plotted using DIAL.  rigure K.5 is the 
full quarter panel with strain (Eps-Y) parallel to the applied 
load. Figures K.6 through K.8 (Eps-Y, Eps-X and Eps-XY) show tne 


Strains in detail close to the cutout. 
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PRNEL RH22 STRIP REINFORCEMENT 
РАМЕ MES ۶۳۶ 





REINFORCEMENT (2 plies/facesheet) 


Figure K.l1 Panel RH22: DIAL Finite Element Mesh. 
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Strain Comparison Around the Cutout. 


Panel RH22 


Figure K.2 
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TABLE XL 


PANEL HH22: LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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Strain Comparison Alonq the X Axis. 


Panel RH22 


Figure K.3 
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TABLE XLI 


PANEL RH22, LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (-10,000 PSI). 
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Microstrain vs. Compressive Stress. 
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Figure K.4 
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TABLE XLII 


SELECTED STRAIN GAGE VALUES DURING LOAD. 


PANEL RH22: 
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Panel RH22 


Figure K.5 
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Eps-Y FEA Contours Near the Cutout. 
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Panel RH22 


Figure K.6 
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Panel RH22 


Figure K.7 
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Panel RH22 


Figure K.8 


APPENDIX L 


PANEL RH31: ANALYTICAL AND EXPERIMENTAL DATA 


Panel RH31 was reinforced with one co-cured ply of G/Ep in the 
shape of two strips on either side of the cutout on the outside of 
each facesheet offset 0.50 inch From the edge of the cutout. The 


reinforcement had the following dimensions: 


Shape: Заир 
Lengtn: 4.70 in 
Width: 1.00 in 


Tnickness (each): С 220522 
Area (each face): 9.40 in 
Total Volume: 05268 in? 


Net Cross Section: 0.056 in? 


The panel failed by facesheet separation and buckling (Туре-2) 


at an applied normal stress of -21,500 psi (O ). Based strictly 


a 
оп the failure of the unreinforced oanel and the computed stress 


concentration factor of 2.82, the predicted failure was at O n" Е 
-32,400 psi. Tne panel failed at 33.0% less than the predicted 
ultimate load. 

The finite element model (mesh) is shown in Figure L.l. The 
area of the strip reinforcement is outlined by the heavy Lines 
0.50 inch to the right of the cutout's edge. 

Figure L.2 compares the three (finite element) computed 


strains around the cutout between the unreinforced panel (PO) 


and RH3l. These computed strain values are listed in Table XLIII. 


259 


There was only a small decrease in the computed strain near the 
cutout due to the reinforcement. 

Figure L.3 compares the computed (solid and dashed lines) and 
experimental (triangles) strains in the X and Y (poisson) direc- 
tions in the panel and shows almost perfect correlation between 
analytical and experimentally measured strain at -10,000 psi 
applied normal stress. There was virtually no strain variation 
between the left and right side of the hole. The edge of the 
reinforcement can not be seen in the figure. The effect of rein- 
forcement is a relatively small decrease in strain from the edge 
of the cutout out to about x = 3.0". Table XLIV gives the 
computed values of the strains along the X axis. 

Figure L.4 shows graphically the stress-strain state during 
the load sequence from om = 0 to -21 ksi. Strain gage values are 
Given an Table XLV: Бо te С. = -12 ksl all gages indicated a 
fairly normal stress-strain relation. At that point up to -18 ksi 
gage #2 demonstrated virtually no strain increase. At -l8 ksi 
gage #1 (x = +0.572") suddenly indicated a strain decrease from 
3900 to o200Me and then an l1004€ increase at -19 ksi. From 
there it remained steady at about 7250 ЦЕ ко failure at Qa = 
-21,500 ksi. Above -18 ksi gage #2 showed a steady decrease in 
strain which most probably indicated a separation of the facesheet 
from the core directly under the gaye. As the load increased tne 
Facesheet ouckled inore and the indicated strain decreased. Тліз 
oanel shows how useful it would have been to instrument both 
facesneets of the panel to measure load transfer between them. 

Figures L.5 tnrough L.8 show the strain contours ша = -19,000 
psi computed and plotted using DIAL. Figure L.5 is the full 


quarter panel with strain (Eps-Y) parallel to the applied load. 
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Figures L.6 through L.8 (Eps-Y, Eps-X and Eps-XY) show the strains 


in detail close to the cutout. 
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C — C — REINFORCEMENT (1 ply/facesheet) 


Figure L.l Panel RH31: DIAL Finite Element Mesh. 
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PANEL RH31: LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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SELECTED STRAIN GAGE VALUES DURING LOAD. 
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Panel RH31: Eps-Y FEA Contours. 


Figure L.5 
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Panel RH31: Eps-Y FEA Contours Near the Cutout. 


Figure L.6 
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Figure L.7 
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Eps-XY FEA Contours Near the Cutout. 


Panel RH31 
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APPENDIX M 


PANEL RH42: ANALYTICAL AND EXPERIMENTAL DATA 


Panel RH42 was reinforced with two co-cured plies of G/Ep in 
the snape of two strips on either side of the cutout on the 
outside of each facesneet offset 0.50 inch from the edge of the 


cutout. The reinforcement had the following dimensions: 


Snape: Strip 

Length: 3.14 in 

Width (each): LOOM in 
Thickness (each): Qu MEE? vo ly) 
Area (each face): 5.280 in^ 

Total Volume: Бр е 


СЕЕ ЕЕЕ ота ша? 


The panel failed at the hole edge (Type-1) at an apolied 
normal stress of -36,990 psi E Based strictly on the failure 
O£ tne unreinforced panel and the comouted stress concentration 
factor of 2.64 the predicted failure was 0. = -34,050 psi. The 
28121 failed at 106.83 of the predicted ultimate stress. It 
sustained the highest load of any test specimen. The reinforce- 
ment increased the panel's weight little more Епап 3.0% and 
increased the failure strength by 21% over the unreinforced panel. 
It was one of the panels that led to tne conclusion that several 
layers of reinforcement close to the cutout are more effective 


than spreading it out more thinly over a larger area. 


Tne finite element model (mesh) is shown in Figure М.1. The 
area of the strip reinforcement is outlined by the heavy lines 
beginning 0.5" to the right of the cutout's edge. 

Figure M.2 compares the three (finite element) computed 
strains around the cutout between the unreinforced panel (PO) 
and RH42. These computed strain values are listed in Table XLVI. 
There is a relatively small decrease in the strain due to the 
strip reinforcement compared with the equivalent amount of rein- 
forcement concentrated next to the cutout. 

Figure M.3 coimoares the computed (solid and dashed lines) and 
experimental (triangles) strains in the Y and X (poisson) direc- 
tions in the panel and shows an excellent correlation between 
analytical and experimental strain at O, = -10,000 psi. There was 
some strain variation between the left and right sides of the 
hole. The exact edge of the reinforcement can not pe seen in the 
figure. The effects of reinforcement is a significant strain 
Aecrease (compared to the unreinforced panel, PO@Y) under the 
reinforcement, a small decrease near the hole anda slight 
increase in the far field (wnere x 2.50"). Table XLVII gives the 
computed values of the strains along the X axis. 

Figure M.4 shows the stress-strain state during the load 
sequence from O, = 0 to -36 ksi. Experimentally measured strain 
data are given in Table XLVIII. Up tO about G = -24 £51 almi 
gages indicated a normal stress-strain relation. At that load 


jage #2 (х = -0.571") demonstrated the "stair-steo" phenomena. At 
-25 ksi gage #1 (х = +0.569") indicates wnat appears to be a 
softening or loss of stiffness--the strain rate drastically 


increased. Gage i2 seems to reflect the same behavior at -4 ksi 


nigher stress. Gage #3 (x = +0.701") appears to pick up the load 
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when gage £l shows what appears to be local ouckling. Note that 
gages #8 and #13 (x = -1.512" and -2.460") reflect none of what is 
Securring next to the cutout. 

The reaction of this panel may help explain much of what 
occurs in the boundary region around the cutout on the other 
panels. At high levels of strain (8,000 to 10,0004€) next to the 
cutout's edge, local delamination, buckling and fiber failure 
forces the transfer of the load path laterally away from the edge 
to the still intact and stiffer fibers and matrix farther from the 
cutout. 

Figures M.5 through M.8 show the strain contours at Gu = 
-10,000 psi computed and plotted using DIAL. Figure M.5 is the 
full quarter panel with strain (EÉps-Y) parallel to the applied 
load. Figures M.6 through M.8 (Eps-Y, Eps-X and Eps-XY) show the 


EE usn detail close to the cutout. 
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Figure M.l Panel RH42: DIAL Finite Element Mesh. 
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Strain Comparison Around the Cutout. 


Panel RH42 
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Figure M. 
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TABLE XLVI 


LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 


PANEL RH42 
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PANEL RH42: LEFEA STRAIN DISTRIBUTION ALONG THE X AXIS (710,000 PSI). 
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Figure M.6 Panel RH42: Eps-Y FEA Contours Near the Cutout. | 
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Figure M.8 


APPENDIX N 


PANEL RH51: ANALYTICAL AND EXPERIMENTAL DATA 


Panel RH51 was reinforced with one co-cured ply of G/Ep in the 
shape of two strips on either side of the cutout on the outside of 
each facesheet offset 0.50 inch from the edge of the cutout. The 


reinforcement had the following dimensions: 


Shape: Sterio 

Length: 7.86 in 

Width (each): “ ۹ 71 
Thickness (each): 8 YQ оу) 


Area (each face): 15:720 in? 
Total Volume: 0.440 in? 


Net Cross Section: 0.056 in? 


The panel failed at the hole edge (Type-1) at an apolied 
normal stress (9) ОЕ -31,030 рз1. Based strictly on the failure 
ОҒ the unreinforced panel and the computed stress concentration 
factor of 2.83, the predicted failure was о" = -32,300 psi. The 
panel failed within 2.1% of the predicted ultimate load. 

The finite element model (mesh) is shown in Figure N.l. The 
area of the strip reinforcement is outlined by the neavy lines 
beginning 0.50 inch to the right of the cutout's edge. 

Figure W.2 compares the three (finite element) computed 
strains around the cutout between the unreinforced panel (POG) 


and RH51. These computed strain values are listed in Table XLIX. 
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There is only a small apparent decrease in.the strain around the 
cutout due to the reinforcement. 

Figure N.3 compares the computed (solid and dashed lines) and 
experimental (triangles) strains in the Y and X (poisson) direc- 
tions in the panel and shows a poorer than usual correlation 
between analytical and experimental strain at om = -10,000 psi. 
There was no apparent strain variation between the left and rignt 
side of the hole, but the finite element model predicted a higher 


level of strain at Or = -10 ksi. From the appearance of the 
panel during the load sequence, this can not be explained. The 
finite element model was rerun to verify the results and the data 
is consistent with the other models. The effect of reinforcement 
is a very slight strain decrease (compared to the unreinforced 
panel, PUO) under the reinforcement and a small decrease near the 
cutout. Table L gives the LZFEA computed values of tne strains 
along the X axis. 

Figure N.4 snows the stress-strain relation during the load 
sequence froin С. = О to -3l ksi. Experimentally measured strain 
values are given in Table LI. Botn gaqes #1 апа #2 (x = +0.509" 
and -0.571") show much higher strain rate than equivalent gages on 
other RH panels. Gage #3 at first parallels tne strain rate of #1 
and #2, then seems to indicate a load transfer away at -3 ksi and 
then again oicks up the load at -9 ksi. Gages #1 and %2 show 
somewhat the same ohenomena described in Appendix M: significant 
Duckling and fiber failure close to the cutout and a trans r 
the Load path away From the cutout's edge. There may have also 
been a transfer of load from one side of the cutout to tne other 


and to the opposing facesheet. 
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Figures N.5 through N.8 show the strain contours at J 


n 

Figure N.5 is the 
full quarter panel with strain (Eps-Y) parallel to the applied 
load. 


-10,000 psi computed and plotted using DIAL. 


Figures N.6 through N.8 (Eps-Y, Eps-X and Eps-XY) show the 


strains in detail close to the cutout. 
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TABLE XLIX 


LEFEA STRAIN DISTRIBUTION AROUND THE CUTOUT (-10,000 PSI). 
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Figure N.5 
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Figure N.7 
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Figure N.8 
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APPENDIX O 


PANEL RR22:  THREE-DIMENSIONAL LINEAR FINITE ELEMENT ANALYSIS 


Panel RR22 was reinforced with two co-cured round plies of 
G/Ep around the cutout on the outside of each facesheet. The 
reinforcement configuration should have been among the most 
efficient, concentrating the maximum amount of reinforcement close 
to the cutout. Figure 5.3 and Table VI show that the round, 200$ 
reinforcement produced about a 22.5$ reduction in maximum strain 
(eps-y) parallel to the applied load. The round, 400$ reinforce- 
ment with twice the volume of additional weight provided only 3.13 
additional strain reduction. It was therefore more than a little 
disconcerting when the most promising panel failed at O, = 21,050 
psi, only 55% of the predicted load. Table IX gives the predicted 
failure (based on the actual failure of the unreinforced panels 
and the LEFEA computed SCF) at 38,250 psi. 

When trying to explain the failure, it was postulated that the 
facesheet layup [03,445,90] may have caused out-of-plane stresses 
(* 0,) sufficient to cause the facesheet to separate from the 
core. This, of course, would have invalidated the entire thesis 
that local reinforcement around a cutout could be a significant 
design benefit. The two-dimensional LEFEA (see section III C.1) 
used in the computational analysis was not able to give stress or 
strain in the Z direction. 

A three-dimensional analysis was undertaken. Figure O.l shows 
the three-dimensional mesh. In order to conserve computer time 


and provide an accurate solution, the quarter panel was modeled 


SOIL 


only from the midplane (z - 0.0). Modeling only half the core and 
and one facesheet did not affect the accuracy of the solution. In 
order to approximate the strain closer to the predicted failure, 
the model was subjected to an equivalent applied load of 30 ksi 
rather than the 10 Ksi used on the 2-D models. The analysis 

was linear and did not take into account the very probable matrix 
cracking and non-linear behavior at high strain ( 10,000 Де). 


Table LII summarizes the results of the analysis. 


Direction: Ү Х 2 XY YZ ZX 
Maximum 1140 4690 22 2420 293 2100 
Minimum -10100 -2720 -11 -9300 -1060 -4220 





Figure 0.2 gives the strain parallel to the applied load. The 
maximum predicted stain was 10,1004€& at О. = 30 ksi. This 15 
exactly three tlmes the maximum strain computed in the 2-D model 
described in Appendix C (see Figure C.5). The exact analytical 
correspondence of the 2- and 3-D FEA helps to validate it. Figure 
O.3 shows the e strain near the cutout. Figure 0.4 shows €. 
next to the cutout. This corresponds with Figure C.7. 

The strain in the Z direction is shown in Figures O.5 and O.6. 
The maximum was 924€, the minimum -llwé. The stress at the 
interface of the facesheet and honeycomb core is shown in Figures 
О.7 апа 0.8. І is obvious that the out-of-plane stress at the 


interface is virtually nil (less than +5 psi) and that premature 
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fallure was not due to the layup or the reinforcement 
configuration. 

The shear strains E € uz 
O.10 and O.ll respectively. The three-dimensional analysis 


and €,. are shown in Figures 0.9, 


reversed the sign on the shear strain from the two-dimensional. 
Comparing the results of the 2- and 3-D analyses, the maximum and 
minimum Є ху: 2420 апа -9300 u€ іп the 3-D (Figure O.9) are almost 
exactly 3 times the 2-D: 826 and -3llOg€ (Figure C.8). 
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Panel RR22: 3-D DIAL Finite Element Mesh. 


Figure О.1 
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Panel RR22: 3-D Eps-Y FEA Contours Near the Cutout. 


Figure O.3 
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Figure O.5 
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APPENDIX P 


FORTRAN PROGRAM “КВЗЕМ“ 


This program written in FORTRAN was developed by S.P. Garbo 
and J.M. Ogonowski of McDonnell Aircraft Company, McDonnell- 
Douglas Aircraft Corporation, PO Box 516, St. Louis, MO 63166. It 
was published by the Air Force Flight Dynamics Laboratory, Wright 
Aeronautical Laboratories, Wright-Patterson Air Force Base, Ohio 
45433 as report AFWAL-TR-81-3041, Volumes 1-3. 

The program was modified by the author to run on the IBM 370. 
The code was renumbered, the input method and output format was 


altered for easier input. 
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90 CONTINUE 
АІ РНА= - АІРНА*РІ/ 180.0 
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DO 160 JJ=1,NUMSTFE 
DO 150 NN NUMPT 


JJJ2JJ-1 
D MU rA EN 105 0 
RADIUSSJJJ*STEINK*DIA/2.0 


С 
<  САБСОБАТЕ x AND Y COORDINATES OF ECINTS ARCUND LOADED HOLE = 


X=RADI US*COS (THETA*+ALPHA 
=RADIUS*SIN (THETA+ALPHA 
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alax ERI 
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110 CONTINUE 
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120 CONTINUE 
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130 CONTINUE 
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1222582 AX ÉSSTY (ATSE] COS (A ALDH 
STRESS (2,d3,NN) =STRX*SLN (ALPH) *SIN(ALPH) +STRY*COS (ALPH) *COS (ALPH) + 
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CALCULATE THE STRESS IN EACH PLY 
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